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ABSTRACT 


The  Point  Sur  Transect  (POST)  was  established  in  1987  by  the  Department 
of  Oceanography  at  the  Naval  Postgraduate  School  to  further  the  present 
understanding  of  long  term  variability  of  current  and  temperature  in  eastern 
boundary  regions.  Of  particular  interest  is  the  temporal  variability  of  poleward 
flows,  their  role  in  gyre-scale  processes,  and  a  more  complete  understanding  of 
their  dynamics.  The  POST  extends  offshore,  normal  to  bottom  topography, 
along  36°  20'N,  to  123°  01.7’W  where  it  doglegs  southwest  along  the  California 
Cooperative  Oceanic  Fisheries  Investigation  (CalCOFI)  line  67.  Station  spacing 
along  the  POST  is  5-10  km  inshore,  increasing  to  20-40  km  offshore.  The 
transect  has  been  occupied  6-8  times  per  year  since  1988,  to  resolve  the  flow  at 
seasonal  and  interannual  time  scales. 

Seven  of  these  cruises  were  selected  for  seasonal  comparisons  of  alongshore 
geostrophic  velocities  and  water  mass  characteristics.  Geostrophic  velocities 
were  referenced  to  the  1000  decibar  surface  based  on  available  PEGASUS 
observations.  Anomalies  of  spiciness  calculated  as  deviations  from  an  average 
offshore  T-S  profile  have  been  used  to  identify  the  location  and  spatial  extent  of 
the  eastern  Pacific  water  masses.  The  California  Undercurrent  (positive 
spiciness  anomaly)  was  a  prominent  feature  in  6  of  the  7  sections  analyzed,  and 
was  very  weak  during  a  period  of  uncommonly  strong  equatorward  wind  stress. 
The  position  of  the  core  varied  from  12  to  42  km  from  shore  while  its  strength 
varied  from  10  to  35  cm  s1,  with  the  maximum  flow  occurring  in  winter.  The 
vertical  extent  of  the  Undercurrent  core  covered  a  region  of  the  continental  slope 
from  70  to  460  m  throughout  these  seven  cruises.  The  nature  of  the  alongshore 
geostrophic  velocities  and  the  location  and  spatial  extent  of  the  undercurrent 
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appear  strongly  related  to  specific  wind  events,  both  local  and  remote.  Remote 
wind  forcing  from  the  south  was  believed  to  cause  anomalous,  strong  poleward 

flow  throughout  the  entire  water  column  during  a  period  of  local  equatorward 

* 

wind  stress,  while  an  intrusion  of  warm,  salty  water  was  believed  to  cause  a  deep 
penetration  (to  700  m)  of  the  California  Current  in  winter. 

These  observations  revealed  primarily  interannual  rather  than  seasonal 
variability.  This  is  not  surprising  since  only  seven  cruises  were  selected  for 
study  whereas  earlier  studies  utilized  many  years  of  data  collected  along  the 
CalCOFI  sampling  grid  to  determine  the  seasonal  means.  This  study  excelled 
over  the  continental  shelf  and  slope  where  the  station  spacing  of  the  POST  is 
considerably  closer  than  the  CalCOFI  scheme,  which  allowed  for  the  study  of 
narrow  coastal  jets  not  well  resolved  by  the  CalCOFI  grid. 
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I.  INTRODUCTION 


Large  scale  atmospheric  forcing  in  the  eastern  Pacific'Ocean  consists  of  the 
North  Pacific  (sub-tropical)  high,  the  Aleutian  low,  and  in  the  summer  the 
thermal  low  over  the  western  United  States.  The  North  Pacific  high  is  most 
intense  during  the  summer  months  while  the  Aleutian  low  is  most  intense  during 
the  winter  months.  The  high  migrates  annually  from  a  maximum  southern 
position  at  28°N,  130°W  in  February  to  a  maximum  northern  position  at  38°N, 
150°W  in  August  (Huyer  1983).  The  U.S.  thermal  low  is  centered  near  35°N 
and  enhances  the  equatorward  wind  stress  over  the  coastal  waters  off  northern 
California  (Reinecker  and  Ehret  1988).  A  region  of  positive  wind  stress  curl 
exists  near  the  coast  throughout  the  year,  being  well  developed  from  May  to 
September  and  having  greater  spatial  variability  during  winter  (Nelson  1977). 

This  large-scale  atmospheric  forcing  creates  the  anticyclonic  North  Pacific 
gyre.  Its  northern  side  is  comprised  of  the  West  Wind  Drift  and  the  North 
Pacific  Current  which  flow  easterly.  The  eastern  limb  of  this  gyre  is  the 
California  Current  (CC).  Offshore  the  CC  is  a  surface  current  (0-300  m  deep) 
carrying  water  equatorward  throughout  the  year  along  the  west  coast  of  North 
America  (Lynn  and  Simpson  1987).  Near  20°N  it  turns  westward  as  part  of  the 
North  Equatorial  Current.  The  average  speed  of  the  CC  off  the  coast  of 
California  is  generally  less  than  25  cm  s1  (Reid  and  Schwartzlose  1962). 

Within  150  km  of  the  coast  there  is  a  fall-winter  reversal  of  the  surface  flow 
known  as  the  California  Countercurrent  (CCC,  Simpson  et  al.  1986)  or  Inshore 
Countercurrent  (IC,  Lynn  and  Simpson  1987).  This  poleward  flow  from 
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October  to  March  is  generally  referred  to  as  the  Davidson  Current  (DC)  north  of 
Point  Conception.  Hereinafter  it  will  be  referred  to  as  the  DC. 

The  California  Undercurrent  (CUC)  flows  poleward  throughout  the  year.  It 
has  been  observed  from  off  Baja  California  (Wooster  and  Jones  1970)  to  as  far 
north  as  Oregon  (Halpem  et  al.  1978).  The  continuity  of  the  CUC  has  not  been 
observed.  It  has  its  origin  in  the  eastern  equator!.-  Pacific  and  is  centered 
primarily  over  the  continental  slope.  The  location,  strength  and  core  depth  as 
determined  from  moored  current  meters  and  inferred  from  hydrographic 
measurements  show  considerable  seasonal  variability  and  can  be  related  to  the 
seasonal  variability  in  wind  stress  and  wind  stress  curl  (Hickey  1979). 

Collectively  these  currents  comprise  what  is  known  as  the  California  Current 
System  (CCS).  The  water  properties  which  make  up  this  system  are  determined 
by  four  water  masses,  each  of  which  can  be  defined  by  its  temperature  (T), 
salinity  (S),  dissolved  oxygen  (DO)  and  nutrient  content  as  it  enters  into  the  CCS. 
Pacific  Subarctic  water  is  formed  at  the  Subarctic  Convergence  through  mixing 
of  the  warmer,  saline  waters  of  the  Kuroshio  Extension  and  the  cooler,  fresher 
waters  of  the  Oyashio  Current  and  enters  the  CC  near  48°N  (Pickard  and  Emery 
1982).  It  is  characterized  by  relatively  low  temperature,  low  salinity,  high 
dissolved  oxygen  and  high  nutrients  (Reid  et  al.  1958).  Equatorial  Pacific  water 
forms  in  the  eastern  tropical  Pacific  and  is  characterized  by  relatively  high 
temperature,  high  salinity,  low  dissolved  oxygen  and  high  nutrients.  It  enters  the 
CCS  from  the  south  and  is  carried  northward  by  the  CUC.  Eastern  North 
Pacific  Central  enters  the  CCS  from  the  west  and  is  characterized  by  relatively 
high  temperature,  high  salinity,  low  dissolved  oxygen  and  low  nutrients  (Reid  et 
al.  1958).  Upwelled  water  within  50  km  of  the  coast  is  identified  by  relatively 
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cold  temperature,  high  salinity,  high  nutrients  and  low  dissolved  oxygen 
(Sverdrup  1938;  Reid  et  al.  1958).  These  water  properties  appear  in  varying 
proportions  throughout  the  CCS  due  to  mixing  which  occurs  as  the  CC  flows 
southward  and  the  CUC  northward. 

Wickham  (1975)  used  5  km  station  spacing  off  Point  Sur  and  Monterey  Bay 
to  define  very  narrow  streams  of  anomalously  high  temperature  and  salinity, 
which  he  termed  southern  water.  Vertical  distributions  of  this  southern  water 
were  extremely  complex.  He  noted  two  regions  of  warmer  southern  water  in 
bands  10-20  km  in  width,  which  appeared  as  intrusions  from  the  south.  Vertical 
cross-sections  indicted  that  these  intrusions  are  generally  found  between  2o0o00 
m  depths.  Alongshore  geostrophic  velocities  along  36°  20’N  exhibited  a  banded 
structure.  A  narrow  band  of  poleward  flow  was  found  between  50  m  and  200  m 
near  the  shelf  edge.  Offshore  he  observed  alternating  bands  (5  km  width)  of 
poleward  and  equatorward  flow,  with  maximum  speeds  of  +40  cm  s'1  and  -80  cm 
S'1,  respectively.  Current  patterns  exhibited  by  both  drogue  measurements  and 
geostrophy  were  in  agreement  with  the  banded  nature  of  southern  water  flowing 
poleward. 

Chelton  (1984)  and  Lynn  and  Simpson  (1987)  have  utilized  the  California 
Cooperative  Oceanic  Fisheries  Investigations  (CalCOFI)  data  sets  to  examine  the 
seasonal  variability  of  alongshore  geostrophic  currents  and  physical  character¬ 
istics,  respectively.  This  large  scale  hydrographic  sampling  grid  was  initiated  in 
1949  and  since  1950  stations  have  been  occupied  between  4  and  12  times  per  year 
(Chelton  1984).  The  sampling  grid  consists  of  a  series  of  parallel  lines  which  are 
oriented  normal  to  the  coast.  Standard  station  and  line  spacing  are  both  74  km, 
with  station  spacing  decreasing  to  40  km  near  the  coast  (Lynn  and  Simpson 


1987).  While  this  data  set  provides  excellent  temporal  coverage,  the  coarse 
station  spacing  and  limited  vertical  extent  (upper  500  m)  do  not  allow  for 
detailed  analysis  of  the  smaller  scale  phenomenon  which  occur  in  coastal  regions. 
The  seasonal  average  geostrophic  flow  in  the  upper  100  m  relative  to  500  m  off 
Point  Sur  and  Point  Conception  is  equatorward  flow  from  February  to 
September  and  poleward  from  October  to  January.  The  deeper  flow  (below  100 
m)  is  different  along  these  sections  with  poleward  flow  throughout  the  year  at 
Point  Conception,  while  off  Point  Sur  it  has  a  poleward  maximum  in  December 
and  weak  equatorward  flow  from  March  to  May.  The  surface  flow  throughout 
this  region  was  found  to  lead  the  annual  wind  forcing,  computed  from  the  spatial 
averages  of  Nelson  (1977),  by  about  one  month  while  the  deep  poleward  flow 
was  found  to  lag  the  local  poleward  barotropic  pressure  gradient  by  about  two 
months  (Chelton  1984).  Chelton  (1984)  estimates  the  total  alongshore  pressure 


gradient  (py  =  barotropic  component  +  baroclinic  component)  as 

f0 


dy 


Pog^  +  g 

dy 


^dz 

dy 


where  the  p  is  pressure,  g  the  gravitational  acceleration,  h  the  sea  surface 
elevation,  D  is  depth,  p  the  water  density,  and  po  the  water  density  at  the  surface. 

Analysis  of  dynamic  height  Fields  have  led  to  the  definition  of  three  domains 
within  the  CCS:  oceanic,  coastal,  and  an  intervening  transition  zone.  The 
transition  zone  is  coincident  with  the  core  of  the  CC  and  is  hypothesized  to  result 
from  recurrent  eddies  and  energetic  meanders.  Seasonal  variability  in  the  fields 
of  temperature,  salinity,  ct  (defined  as  at  =  ps>t>o  -  1000,  where  p  has  units  of  kg 
nr3),  and  oxygen  is  related  to  variations  in  the  CC,  DC,  and  CUC  through 
vertical  adjustments  in  the  density  field  and  through  changes  in  transport  (Lynn 
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and  Simpson  1987).  This  transition  zone  differs  from  that  described  in  the 
Coastal  Transition  Zone  (CTZ)  program  where  it  is  defined  as  the  area 
encompassing  the  offshore  (greater  than  50  km)  region  of  cold  "filaments”  often 
seen  in  satellite  imagery  (Brink  and  Hartwig  1985). 

The  Central  California  Coastal  Circulation  Study  (CCCCS)  (February  1984 
and  July  1984)  provided  higher  resolution  (than  CalCOFI)  CTD  and  current 
meter  coverage  from  San  Francisco  to  Point  Conception  and  was  aimed  at  a 
detailed  description  of  the  coastal  circulation  on  the  continental  shelf  and  upper 
slope  (Chelton  et  al.  1988).  During  the  first  half  of  1984  the  mean  flow  between 
Point  Conception  and  Point  Sur  was  poleward  and  in  opposition  to  the 
equatorward  wind  stress.  A  three  week  period  of  calm  winds  during  July  1984 
resulted  in  a  100  km  wide  surface  poleward  flow  (observed  in  current  meter 
measurements,  satellite  sea  surface  temperautre  imagery,  buoy  drift  trajectories, 
and  inferred  from  the  hydrographic  data)  which  extended  300  km  off  the  central 
California  coast.  In  July  1981  a  similar  but  weaker  occurrence  produced  a 
poleward  surface  flow  which  extended  to  150  km  off  the  coast.  Poleward  flow 
at  the  surface  over  the  continental  shelf  is  normally  observed  during  spring  and 
summer,  however,  poleward  flow  at  the  surface  over  the  continental  slope  was 
previously  not  observed  after  February  (Chelton  et  al.  1988). 

In  1987,  the  Point  Sur  Transect  (hereinafter  referred  to  as  the  POST)  was 
established  by  the  Department  of  Oceanography  at  the  Naval  Postgraduate  School 
(NPS)  to  further  the  present  understanding  of  long  term  variability  in  eastern 
boundary  regions.  Of  particular  interest  is  the  time  variability  of  poleward 
flows,  their  role  in  gyre-scale  processes,  and  a  more  complete  understanding  of 
the  dynamics  involved. 


The  purpose  of  this  study  is  to  describe  the  seasonal  variability  of  the 
alongshore  geostrophic  currents  and  water  mass  structure  off  Point  Sur  using 
hydrographic  data  four  times  per  year  collected  from  spring  1988  through 
winter  1989  along  the  POST.  Data  collection  and  processing  will  be  discussed  in 
Chapters  II  and  III.  Analysis  of  water  mass  characteristics  and  alongshore 
geostrophic  velocities  are  discussed  in  Chapter  IV  and  comparisons  with  earlier 
studies  conducted  off  Point  Sur  are  made  in  Chapter  V.  A  discussion  and 
conclusion  follows. 
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II.  DATA  COLLECTION  AND  PROCESSING 


The  POST,  as  established  in  1987,  extends  offshore,  normal  to  bottom 
topography,  along  36°  20'N,  to  123°  01.7'W  where  it  meets  and  thereafter 
coincides  with  the  California  Cooperative  Oceanic  Fisheries  Investigation 
(CalCOFI)  line  67  (the  Monterey  Bay  line).  It  follows  line  67  to  allow  for 
comparison  with  previous  studies  which  have  utilized  the  CalCOFI  data  set.  As 
of  January  1990  there  have  been  16  cruises  (Table  1)  along  the  POST  using 
either  one  of  two  sampling  schemes.  The  first,  used  for  the  NPS  student  cruises, 
consists  of  19  stations  and  extends  to  approximately  370  km  offshore  (Figure 
la).  The  second,  used  for  CUC  cruises,  consists  of  22  stations  and  is 
approximately  215  km  in  length,  terminating  in  the  vicinity  of  station  15  of  the 
NPS  student  cruise  sampling  scheme  (Figure  lb).  The  former  provides  greater 
offshore  coverage  while  the  latter  provides  a  more  extensive  coverage  of  the 
inshore  end  of  the  transect  for  a  more  detailed  study  of  the  CUC.  Both  of  these 
schemes  provide  higher  resolution  coverage  than  the  CalCOFI  sampling  scheme 
in  this  region  (Figure  2).  Of  the  3  6  cruises  along  the  POST,  a  total  of  7, 
indicated  by  asterisks  in  Table  1,  were  selected  for  detailed  analysis.  These 
cruises  provide  a  suitable  seasonal  description  of  the  POST.  When  cruises  were 
within  a  month  of  one  another  only  one  was  selected. 

A.  DATA  COLLECTION 

Data  collected  along  the  POST  consists  primarily  of  three  types: 
hydrographic  (CTD),  acoustic  Doppler  current  profiler  (ADCP),  and  more 
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I  or  I 


Figure  1.  Sampling  schemes  along  the  Point  Sur  Transect: 

(a)  NFS  student  cruise  sampling  scheme;  (b)  California 
undercurrent  cruise  sampling  scheme. 
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125°  120° 


Figure  2.  CalCOFI  sampling  scheme  along  the  central  California 
coast:  Dots  represent  hydrographic  stations  along  cardinal 
lines  which  lie  in  water  deeper  than  5(X)m  and  were  occupied 
r.iore  than  40  times  between  1950  and  1979.  (Source:  Chelton  1984) 
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Table  1.  POINT  SUR  TRANSECT  CRUISE  PERIODS  AND  DATA 

TYPES:  Includes  all  cruises  through  January  1990.  Cruises  used  in  this 
study  are  indicated  by  asterisks. 


Cruise 

Dates 

Vessel 

Data  Type 

STNOV  1987 

11/4  -  11/11 

R/V  Point  Sur 

CTD,  ADCP 

CUC-April  1988 

4/15-  5/1 

USNS  DeSteiguer 

CTD,  ADCP, 
PEGASUS 

STM  AY  1988  * 

5/4  -  5/11 

R/V  Point  Sur 

CTD,  ADCP 

CUC- August  1988  * 

8/3  -  8/7 

R/V  Point  Sur 

CTD,  ADCP, 
PEGASUS 

CUC-September  1988 

9/22  -  9/27 

R/V  Point  Sur 

CTD,  ADCP, 
PEGASUS 

STNOV  1988 

11/1  -  11/8 

R/V  Point  Sur 

CTD,  ADCP 

CUC-November  1988  * 

11/14  -  11/19 

R/V  Point  Sur 

CTD,  ADCP, 
PEGASUS 

CUC-February  1989  * 

2/3  -  in 

R/V  Point  Sur 

CTD,  ADCP, 
PEGASUS 

CUC-March  1989 

3/24  -  3/30 

R/V  Point  Sur 

CTD,  ADCP, 
PEGASUS 

CUC- May  1989 

5/10  -  5/26 

USNS  DeSteiguer 

CTD,  ADCP, 
PEGASUS 

STM  AY  1989  * 

5/4  -  5/8 

R/V  Point  Sur 

CTD,  ADCP 

CUC-July  1989  * 

7/28  -  8/3 

R/V  Point  Sur 

CTD,  ADCP, 
PEGASUS 

CUC-September  1989 

9/25  -  9/30 

R/V  Point  Sur 

CTD.  ADCP, 
PEGASUS 

STNOV  1989 

11/1  -  11/8 

R/V  Point  Sur 

CTD.  ADCP 

CUC-November  1989  * 

11/15  -  11/22 

R/V  Point  Sur 

CTD.  ADCP, 
PEGASUS 

CUC-January  1990 

1/17  -  1/24 

R/V  Point  Sur 

CTD,  ADCP, 
PEGASUS 

10 


recently  data  collected  using  the  PEGASUS,  a  free-falling  acoustically  tracked 
dropsonde  that  measures  temperature,  salinity,  velocity,  and  pressure.  To 
provide  additional  temperature  information,  expendable  bathytnermographs 
(XBT's)  were  also  utilized  on  some  student  cruises.  A  listing  of  the  data 
collected  on  each  individual  cruise,  including  the  vessel  used,  can  be  found  in 
Table  1.  Data  analysis  in  this  thesis  has  been  limited  to  hydrographic 
information  from  CTD  casts. 

Hydrographic  data  was  collected  using  Neil  Brown  Mark  III-B  CTDs.  The 
only  distinction  between  these  two  instruments  is  the  physical  size  of  their 
pressure  casings.  Because  of  this,  they  have  been  commonly  referred  to  as  either 
the  "large"  or  "small  "  Neil  Brown  CTD,  a  convention  which  will  be  used 
throughout  this  thesis.  The  Neil  Brown  Mark  III-B  CTD  has  a  resolution  of 
±0.001  PSU,  ±0.005°C,  and  0.0015%  of  the  depth  range,  and  is  considered 
accurate  to  within  ±  0.005  PSU,  ±  0.005°C,  and  0.1%  of  the  depth  range  for 
salinity,  temperature,  and  pressure,  respectively.  On  each  cruise,  a  General 
Ocean ics  Rosette  sampler,  equipped  with  twelve  5-liter  bottles,  was  attached  to 
trie  CTD  for  in  situ  water  sampling.  These  samples  were  used  for  post-cruise 
calibration  of  the  hydrographic  data.  Conductivity,  temperature,  and  pressure 
data  were  acquired  on  the  downcast  at  each  station,  while  bottle  samples  were 
collected  on  the  upcast.  Winch  speeds  were  approximately  0.5  m  s  '  for  the 
upper  and  lower  150  m  of  the  cast  and  1.0  m  s_1  in  between. 

Raw  wind  data  from  the  NOAA  weather  buoys  (Figure  3)  located  off 
Monterey  Bay  (B46042)  and  Cape  San  Martin  (B46028)  was  obtained  from  Mr. 
Dave  Husby,  Pacific  Fisheries  Environmental  Group  (PFEG).  Data  gaps  ranging 
from  one  to  several  hours  existed  and  required  visual  interpolation.  With  the 
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exception  of  buoy  46042,  which  was  out  of  service  during  January  and  February 
1989,  complete  wind  records  for  each  cruise  were  obtained. 

B.  DATA  PROCESSING 

The  raw  CTD  data,  which  was  collected  at  0.1  m  intervals,  was 
averaged  in  2  meter  bins  using  programs  written  by  Mr.  Paul  lessen,  NPS 
Oceanography  Dept.  The  data  was  then  calibrated  according  to  the  procedures 
discussed  in  chapter  III.  Additional  quality  control  consisted  of  searching  the 
temperature  and  salinity  information  for  vertical  gradients  in  excess  of  0.2°C  nrr1 
and  1.0  PSU  nr1,  respectively.  When  these  values  were  exceeded,  the  data  was 
visually  inspected,  and  if  it  was  determined  that  the  points  were  in  error,  the  data 
was  linearly  interpolated  (lessen  et  al.  1989).  Errors  of  this  type  are  usually  the 
result  of  the  temporal  mismatch  in  the  response  of  the  temperature  and  salinity 
sensors. 

To  allow  for  the  computation  of  anomalies  of  temperature,  salinity,  and 
spiciness  along  a  constant  pressure  surface,  these  data  sets  were  interpolated  to  2 
dbar  pressure  levels.  This  interpolation  of  the  original  data  sets,  averaged  to  2 
meter  bins,  was  necessary  because  this  averaging  did  not  always  result  in  an 
incremental  value  of  exactly  2  dbar.  Density  anomaly  (y)  was  was  computed 
using  the  algorithm  found  in  Volume  4  of  the  International  Oceanographic 
Tables  (UNESCO,  1987).  Spiciness,  which  is  a  state  variable  most  sensitive  to 
isopycnal  variations,  was  computed  using  the  algorithm  developed  by  Flament 
(1986). 

The  wind  data  from  buoys  46028  and  46042  were  first  transformed  into 
a  u,  v  coordinate  system  and  then  into  an  alongshore  and  across-shore  coordinate 
system.  The  drag  coefficient  was  calculated  (Large  and  Pond  1981)  as 
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Cd  =  1.14  x  10*3,  Ivl^lOms'1 

CD  =  0.49  x  1  O'3  +  (0.065  x  10'3)(|v|),  |v|>  10  ms’1 

The  alongshore  wind  stress  was  the  computed  as 

^AL=  PhCdIvIval 

where  p a is  the  density  of  air,  lv  I  is  the  magnitude  of  the  wind  velocity,  and  val 
is  the  alongshore  component  of  the  wind  velocity. 
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Figure  3.  NOAA  weather  buoys  used  to  calculate  alongshore 

component  of  wind  stress:  Buoy  46042  is  off  Monterey  Bay 
and  buoy  46028  is  off  Cape  San  Martin.  The  100,  200,  5(H),  and 
1  OCX)  m  isobaths  are  also  shown.  (Adapted  from:  Chelton  et  al. 
1988) 
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III.  CALIBRATION  PROCEDURE  FOR  HYDROGRAPHIC 

DATA 

t 

Prior  to  any  type  of  formal  scientific  analysis,  a  data  set  should  be  subject  to 
a  rigorous  calibration  procedure  to  produce  the  most  accurate  results  possible. 
This  section  discusses  the  methodologies  used  in  both  pre-  and  post-  cruise 
conductivity,  temperature  and  pressure  calibrations,  along  with  a  more  detailed 
discussion  of  the  salinity  calibration  procedure  used  in  the  post  processing  of 
individual  cruises  using  bottle  data. 

A.  STANDARD  PRE-  AND  POST-  CRUISE  CALIBRATION 
PROCEDURE 

Prior  to  each  cruise  a  series  of  calibrations  were  performed  to  ensure  the 
accuracy  of  the  CTD's  conductivity,  temperature  and  pressure  sensors.  In  some 
instances,  a  similar  series  of  calibrations  was  performed  after  the  cruise  had  been 
completed.  In  general,  however,  the  pre-cruise  calibration  for  one  cruise  served 
as  the  post-cruise  calibration  for  the  previous  cruise.  All  of  the  calibrations 
discussed  in  this  section  were  carried  out  by  NPS  personnel. 

The  conductivity  calibrations  were  carried  out  using  either  a  Guildline 
Model  8400  Autosal  or  an  AGE  Instruments  Model  2100  Minisal.  Comparisons 
between  the  standard  and  the  CTD  sensor  were  made  at  several  different 
conductivity  levels.  At  each  level,  a  number  of  samples  were  taken  and  averaged 
to  yield  a  single  conductivity  value  for  the  sensor.  Based  upon  the  standard  and 
sensor  readings,  a  regression  analysis  was  used  to  obtain  the  coefficients 
necessary  to  correct  the  sensor  to  the  standard.  In  all  cases  the  best  fit  was  a 
simple  linear  regression.  Since  salinity  and  conductivity  work  equally  well  as  the 
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calibration  parameter,  salinity  was  chosen  for  convenience  based  upon  the 
structure  of  NPS  oceanographic  data  storage. 

The  temperature  calibrations  were  carried  out  using  Seabird  temperature 
sensors  as  the  reference  standard.  Three  different  sensors  (  SBE  #  439,  626,  and 
664)  were  employed,  most  frequently  SBE#  664.  These  sensors  were  calibrated 
by  the  manufacturer  approximately  every  two  years.  In  addition,  there  were 
several  in-house  calibrations  performed  between  the  manufacturers  calibrations. 
In  these  latter  calibrations,  one  Seabird  sensor  was  calibrated  against  another  for 
a  quick  comparative  check  on  the  accuracy.  Although  these  in-house  calibrations 
were  not  as  rigorous  as  those  performed  by  the  manufacturer,  they  still  allow  for 
early  detection  of  large  drifts  in  these  sensors.  For  each  calibration,  an  insulated 
tub  containing  approximately  70  to  80  liters  of  fresh  water  was  used. 
Comparisons  were  made  between  the  reference  standard  and  the  CTD 
temperature  sensor  at  1°C  increments  from  0°C  up  through  30°C.  At  each 
temperature,  a  number  of  values  were  collected  and  averaged  to  yield  a  single 
value.  These  values  were  then  used  to  determine  the  best  fit  to  the  reference 
standard,  with  the  result  again  being  a  simple  linear  regression.  The  values  of 
the  temperature  coefficients  used  for  each  cruise  (Table  2)  are  in  most  cases  the 
pre-cruise  values,  although  in  some  instances,  the  post-cruise  values  were  used. 
Where  the  difference  between  the  pre-  and  post-cruise  temperature  calibrations 
approached  or  exceeded  the  accuracy  of  the  instrument,  the  instrument  and 
reference  temperature  values  from  each  calibration  were  combined  into  one  set 
for  determination  of  the  final  coefficients. 

The  pressure  calibrations  were  performed  with  a  Chandler  Engineering 
deadweight  tester  as  the  reference  standard.  Comparisons  made  between  the 
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CTD's  pressure  sensor  and  the  deadweight  tester  at  various  pressures  yield 
results  that  are  within  the  manufacturers  specifications  of  0.1%  of  full  scale.  In 
some  instances,  a  slope  of  1.0  was  used  while  in  others,  a  value  of  other  than  1.0 
was  used  (Table  2).  In  all  cases,  the  pressure  offset,  which  is  this  observed 
difference  between  the  CTD  pressure  sensor  reading  and  the  actual  pressure  on 
deck  at  the  start  of  a  cast,  was  applied  to  the  data  set  and  serves  as  the  pressure 
intercept. 


Table  2.  TEMPERATURE  AND  PRESSURE  CALIBRATION 
COEFFICIENTS 


Cruise 

Temperature 

Slope 

Temperature 

Intercept 

Pressure 

Slope 

STMAY  1988 

0.998543 

0.047536 

1.000000 

0.998030 

0.039016 

0.999640 

|  November  1988 

0.999769 

0.0097487 

0.999950 

■■■ 

0.999C80 

0.002360 

0.999638 

[  STMAY  1989 

0.999364 

0.003435 

1.000000 

0.999960 

-0.001150 

0.999393 

November  1989 

0.999950 

0.001800 

1.000150 

B.  POST-CRUISE  SALINITY  CALIBRATION  PROCEDURE 

The  previous  section  dealt  strictly  with  the  calibration  procedures  that  are 
conducted  routinely  at  the  beginning  and  end  of  each  cruise.  By  themselves  they 
are  generally  not  sufficient  to  guarantee  an  accurate  data  set.  In  this  section  the 
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manner  in  which  the  bottle  samples  collected  on  each  cruise  were  used  to 
perform  a  calibration  upon  salinity  is  discussed. 

At  the  end  of  each  cruise,  the  bottle  samples  were  analyzed  using  either  the 
Guildline  Model  8400  Autosal,  which  determines  the  conductivity  ratio  between 
the  sample  and  the  reference  standard,  or  the  AGE  Instruments  Model  2100 
Minisal,  which  outputs  the  bottle  salinity.  The  program  used  to  compute  salinity 
from  conductivity,  temperature  and  pressure,  CONDCAL  FORTRAN,  was 
written  by  Mr.  Jim  Stockel  and  utilizes  the  algorithm  of  Lewis  and  Perkin 
(1981).  CONDCAL  corrects  both  the  raw  temperature  and  pressure  for  the 
pre/post  cruise  calibrations,  including  the  pressure  offset  and  produces  a  salinity 
value  corresponding  to  the  input  conductivity  value.  In  addition,  this  program 
will  also  compute  both  the  true  conductivity  and  salinity  based  upon  the 
conductivity  ratio,  provided  that  the  Autosal  was  used.  Once  the  raw  salinity  was 
computed,  it  was  subtracted  from  the  bottle  salinity,  hereinafter  referred  to  as 
true  salinity,  to  obtain  the  salinity  difference,  used  to  determine  how  well  the 
CTD  salinities  have  been  calibrated. 

Before  the  calibration  procedure  could  begin,  it  was  necessary  to  separate 
those  cruises  on  which  the  small  Neil  Brown  CTD  was  used  from  those  on  which 
the  large  Neil  Brown  CTD  was  used.  The  reason  for  this  separation  was  that  the 
small  CTD  conductivity  channel  had  a  peculiar  pressure  dependence,  which  was 
very  noticeable  when  salinity  difference  was  plotted  versus  pressure  (depth). 
The  error  between  the  true  salinity  value  and  that  of  the  instrument  became 
greater  with  increasing  depth.  The  problem  was  traced  to  a  bad  conductivity 
cell,  which  was  replaced  prior  to  the  November  1989  California  Undercurrent 
Cruise.  However,  those  cruises  which  used  the  small  CTD  prior  to  November 
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1989  were  subject  to  this  problem  and  therefore  had  to  be  calibrated  separately. 
Because  of  this  problem,  three  separate  calibration  procedures  were  used;  one 
for  the  large  CTD  and  one  for  the  small  CTD  both  before  ?md  after  repair.  Each 
procedure  will  be  discussed  at  length  in  the  following  sections. 

1 .  Calibration  of  Large  Neil  Brown  CTD 

A  total  of  105  bottles  samples  were  collected  at  83  CTD  stations  for 
STMAY  1988.  These  samples  were  collected  on  two  separate  legs,  the  first  off 
Point  Sur,  and  the  second  in  Monterey  Bay.  With  the  exception  of  one  station  on 
each  leg,  there  was  only  one  bottle  sample  taken  at  the  bottom  of  each  station.  A 
scattergram  of  the  initial  bottle/instrument  differences  (Figure  4)  shows  that 
some  outliers  existed  and  needed  to  be  removed.  The  removal  of  these  outliers  is 
based  upon  the  assumption  that  these  differences  are  randomly  distributed  and 
should  fluctuate  very  little  with  time  and  depth.  These  points  are  due  to  ship 
motion  in  high  gradients,  human  error  in  copying  numbers,  etc.  This  statement 
also  assumes  that  any  temporal  trends  are  small  and  approach  the  noise  level  of 
the  instrument.  For  the  large  CTD,  these  appear  to  be  valid  assumptions. 

The  criteria  chosen  in  most  cases  to  eliminate  outliers  was  a  two 
standard  deviation  limit.  At  most,  two  passes  were  made  through  each  data  set. 
This  choice  was  made  arbitrarily  for  statistical  consistency.  After  each  pass  a 
new  standard  deviation  was  computed  and  used  as  the  criteria  for  the  next  pass. 
Points  that  fell  outside  of  these  limits  were  subject  to  removal  if  they  appeared 
unrealistic  with  respect  to  the  observed  salinity  difference  pattern.  If  points 
were  close  to  the  limit  they  were  retained  for  further  analysis. 
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Figure  4.  Scattergram  of  initial  salinity  difference  for  STMAY 
1988:  This  plot  results  from  subtracting  the  instrument  salinity 
from  the  bottle  salinity.  No  outliers  have  been  removed.  Positive 
values  indicate  the  instrument  is  reading  too  low.  The  mean  (m)  is 
0.828  and  the  standard  deviation  (sd)  is  0.253. 

As  a  result  of  the  first  pass  through  this  cruise  data,  three  points  were 
eliminated  from  further  consideration.  These  are  the  highlighted  points  in 
Figure  4.  The  effect  of  removing  these  points  can  be  seen  in  Figure  5,  which 
also  indicates  the  existence  of  additional  outliers.  Notice  that  the  scale  is  not  as 
coarse  as  that  shown  in  Figure  4  due  to  the  removal  of  these  points.  The 
scattergram  shown  in  Figure  6  results  from  removing  three  more  points,  and  is 
the  data  set  used  to  calibrate  the  salinities  for  this  cruise.  Notice  again  the  change 
to  a  finer  scale  with  the  removal  of  these  points. 
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Figure  5.  Scattergram  of  intermediate  salinity  difference  for 
STMAY  1988:  This  plot  results  from  removing  three 
outliers.  Points  removed  were  the  negative  values  shown  in 
Figure  4.  Positive  values  indicate  the  instrument  is  reading 
too  low.  The  mean  (p)  is  0.868  and  the  standard  deviation  (sd) 
is  0.047. 

After  this  process  was  completed,  a  total  of  99  data  points  remained  for 
use  in  calibrating  the  instrument  relative  to  the  bottle  samples.  The  best  fit  was 
linear  (Figure  7)  with  a  slope  of  1.036714  and  an  intercept  of  -0.347450.  The 
effect  of  this  regression  was  to  increase  the  instrument  salinities,  thereby 
removing  an  almost  constant  offset  of  0.868  PSU.  The  resulting  residual 
statistics  had  a  mean  of  zero  and  a  standard  deviation  of  0.007,  with  a  minimum 
residual  of  -0.032  and  a  maximum  of  0.016.  Here  the  residual  is  computed  as 
the  fitted  instrument  salinity  minus  the  raw  instrument  value.  The  standard 
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error  of  the  residual  mean  was  0.001  and  the  correlation  coefficient  (r^)  was 
0.999588.  This  latter  statistic  is  an  indicator  of  how  much  of  the  total  variance 
can  be  explained  by  a  simple  linear  regression.  In  this  case  it  indicates  a  strong 
linear  relationship. 

The  same  calibration  procedure  was  followed  for  each  of  the  remaining 
cruises  using  the  large  CTD.  The  resulting  calibration  coefficients  can  be  found 
in  Table  3.  In  all  cases,  a  linear  fit  was  used  in  the  final  calibration  of  the 
instrument  values.  This  choice  was  based  upon  the  simple  linear  relation 
exhibited  between  the  instrument  and  true  salinity  values.  Calibration  of  the 
small  CTD  required  considerably  more  effort  and  will  be  discussed  in  the  next 
section. 

2.  Calibration  of  Small  Neil  Brown  CTD  before  repair 

The  small  CTD  was  observed  to  have  an  unusual  pressure  dependence  as 
a  result  of  a  faulty  conductivity  cell  which  produced  larger  instrument  errors  as 
the  ambient  pressure  was  increased.  As  a  result  of  this  problem,  the  calibration 
procedure  required  a  modification  to  remove  this  pressure  effect.  A  second 
order  polynomial  fit  was  used  in  each  case  to  model  the  pressure  dependence. 
Fitted  salinity  differences  were  computed  from  these  polynomials  and  added  to 
the  instrument  salinities  to  remove  the  pressure  effect.  Once  the  pressure  effect 
was  removed,  the  calibration  proceeded  in  a  manner  similar  to  those  for  the 
large  CTD,  and  the  procedure  discussed  has  been  applied  to  all  small  CTD 
cruises  prior  to  November  1989. 

The  discussion  of  this  procedure  has  been  divided  into  two  parts: 
removal  of  the  pressure  effect  and  final  calibration  of  instrument  salinity  values. 
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Figure  6.  Scattergrani  of  final  salinity  difference  for  STMAY  1988: 

This  plot  results  from  removing  three  more  three  outliers.  Points 
removed  were  the  values  below  0.800  in  Figure  5.  Positive 
values  indicate  the  instrument  is  reading  too  low.  The  mean  (|i)  is 
0.875  and  the  standard  deviation  (sd)  is  0.015. 


Figure  7.  Linear  Regression  used  to  ft  STMAY  1988  CTD 
salinities  to  bottle  salinities:  This  plot  results  from 
regressing  STMAY  1988  CTD  salinities  against  corresponding 
bottle  samples.  The  correlation  coefficient  (r^)  is  0.999588,  which 
indicates  that  99.9588%  of  the  variance  can  be  explained  by  a 
simple  linear  regression  model. 
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Table  3.  SALINITY  CALIBRATION  COEFFICIENTS  FOR  THE 
LARGE  NEIL  BROWN  CTD:  The  standard  deviation  listed  in 
this  table  is  of  the  residuals  after  the  calibration  was  applied.  The 
number  of  points  indicates  the  number  of  instrument/bottle  pairs 
used  in  the  analysis. 


Cruise 

Slope 

Intercept 

#  of 
points 

Standard  Deviation 

STMAY  1988 

1.036714 

-0.347450 

99 

±0.007 

1.043071 

-0.568235 

45 

±0.006 

|  STMAY  1989 

1.005162 

-0.167634 

37 

±0.005 

1.002600 

-0.140919 

45 

±0.005 

Removing  this  pressure  effect  proved  to  be  a  very  difficult  task  and 
unfortunately  residual  statistics  were  not  as  good  as  in  the  case  of  the  large  CTD. 
a.  Pressure  effect  and  its  removal 

To  illustrate  this  process,  the  data  from  the  February  1989  cruise 
will  be  discussed  here  at  length.  This  particular  cruise  actually  required  two 
separate  corrections  to  remove  the  pressure  effect  and  was  chosen  for  discussion 
for  this  reason.  The  pressure  effect  (Figure  3)  had  a  paiabolic  shape,  and  in  all 
but  this  particular  case,  a  single  quadratic  polynomial  fit  was  used  to  model  and 
remove  the  pressure  effect.  During  the  course  of  this  cruise  the  shape  of  the 
curve  appeared  to  change.  The  resulting  curve  from  stations  18  through  32  had 
considerably  more  curvature  to  it,  unlike  the  curve  from  earlier  stations.  This  is 
illustrated  in  Figure  8  w^ere  the  lower  set  of  points  belong  to  stations  1  through 
17,  and  the  upper  set  to  stations  18  tl. rough  32.  Attempting  to  use  this  curve  to 
correct  the  data  prove  unsatisfactory  since  predicted  salinities  were  much  too 
high  in  comparison  to  corresponding  bottle  salinity  values.  Using  the  deep  T-S 
curves  as  a  reference,  the  best  results  were  obtained  w'hen  stations  1  through  17 
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were  calibrated  independently,  and  stations  18  through  32  were  calibrated  using 
the  polynomial  computed  from  all  32  stations.  Using  bottle  samples  from  all 
stations  to  calibrate  those  from  stations  18  through  32,  allowed  the  earlier  bottles 

f 

to  influence  the  polynomial  fit,  which  resulted  in  a  shift  of  the  computed  curve 
closer  to  the  characteristic  shape  exhibited  by  this  instrument  on  this  and  other 
cruises  using  the  small  CTD. 

The  procedure  described  earlier  for  removing  the  pressure 
dependence,  utilized  pressure  as  the  input  variable  into  the  quadratic  polynomial 
to  obtain  a  fitted  value  for  the  salinity  difference.  This  fitted  value  was  then 
added  to  the  raw  instrument  salinity  yielding  a  salinity  value  absent  of  the 
pressure  effect.  A  new  salinity  difference  was  computed  by  subtracting  this 
salinity  from  the  true  salinity  and  was  used  in  subsequent  steps  of  the  calibration 
procedure.  The  effectiveness  of  this  procedure  on  both  data  sets  can  be  seen  in 
Figure  9  which  illustrates  the  salinity  difference  before  and  after  removal  of  the 
pressure  effect.  The  results  show  an  improvement,  however,  they  are  not  quite 
as  good  as  those  in  the  case  of  the  large  CTD,  and  in  light  of  the  problems 
discussed  earlier,  they  have  been  taken  to  be  the  best  obtainable. 

b.  Final  adjustment  of  Small  Neil  Brown  CTD  salinities 

The  next  step  was  to  compare  the  pressure  adjusted  salinities  to  the 
true  salinities  and,  if  necessary,  determine  any  further  corrections.  The 
procedure  used  for  this  was  identical  to  that  used  with  the  large  CTD,  namely, 
examine  a  scattergram  of  the  salinity  differences  to  detect  and  remove  any 
outliers,  and  then  use  a  linear  regression  to  fit  the  instrument  salinities  to  the  true 
values. 
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Figure  8.  Pressure  effect  on  small  Neil  Brown  CTD  (February 

1989):  This  plot  illustrates  the  pressure  effect  in  terms  of  salinity 
difference.  Negative  values  for  salinity  difference  indicate  the 
instrument  is  reading  too  high.  Pressure  is  plotted  in  terms  of 
corrected  pressure  based  upon  the  pre-cruise  calibration  and  pressure 
offset. 


Scattergram  of  Salinity  Difference 


o  Sal.  Diff. 

#  New  Sal  Diff 


Figure  9.  Comparison  of  salinity  differences  before  and  after 
removing  pressure  effect  from  CTD  data  (February 
1989):  This  plot  is  results  from  removing  the  pressure 
dependence  from  the  original  salinity  differences.  The  open 
circles  represent  salinity  differences  prior  to  removal,  whereas 
the  blackened  circles  represent  the  new  salinity  differences  after 
the  removal. 


The  final  step  in  the  calibration  of  the  February  1989  cruise  data 
was  to  perform  a  regression  between  instrument  and  true  salinity.  Again,  two 
separate  calibrations  were  used  over  the  same  data  sets.  In  both  cases,  the  chosen 
fit  was  linear  (Figures  10a  and  10b).  Stations  0  through  17  were  calibrated  as 
one  group,  stations  18  through  32  using  all  stations.  Station  0  was  the  first 
station  occupied  on  this  cruise,  however  no  bottle  samples  were  taken. 

The  net  result  of  these  calibrations  were  very  small  shifts  in  the 
instrument  salinities  which  indicates  that  the  pressure  effect  was  the  cause  of 
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most  of  the  observed  difference.  The  choice  of  the  linear  fit  was  again  based 
upon  the  linear  appearance  of  the  data. 

The  calibration  of  the  November  1988  cruise  required  only  a  single 
calibration  to  remove  the  pressure  effect.  Upon  examination  of  the  resulting 
salinity  differences,  it  was  noticed  that  two  distinct  patterns  existed  (see 
Appendix  A).  Stations  1  through  17  were  calibrated  separately  as  one  group, 
while  stations  18  through  32  were  calibrated  as  another.  The  coefficients  used 
for  this  cruise  as  well  as  for  February  1989  can  be  found  in  Tables  4  and  5. 

3.  Calibration  of  Small  Neil  Brown  CTD  after  repair. 

The  cause  of  the  pressure  dependence  associated  with  the  small  CTD 
w'as  a  faulty  conductivity  cell.  Before  the  small  CTD  was  used  on  the  November 
1989  cruise  a  new  conductivity  cell  was  installed  by  NPS  personnel  and  the  unit 
recalibrated.  At  this  time  a  large  portion  of  the  pressure  effect  described  above 
was  removed.  The  calibration  procedure  used  for  the  November  1989  cruise 
data  was  much  simpler  that  earlier  small  CTD  calibrations  as  a  result  of 
instrument  modifications  performed  prior  to  the  cruise.  The  initial  scattergram 
of  salinity  difference  versus  pressure  indicated  that  a  simple  two  step  linear 
pressure  correction  was  required.  In  this  case  one  correction  was  applied  to 
pressures  less  than  or  equal  to  1000  dbar,  while  the  other  w'as  for  pressure 
greater  than  1000  dbar.  To  eliminate  any  discontinuities  in  the  data  both  linear 
equations  were  constrained  to  yield  identical  results  at  the  1000  dbar  level. 
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Figure  10.  Linear  Regression  to  fit  February  1989  CTD  salinities  to 
bottle  salinities:  (a)  represents  the  fit  for  stations  18-32,  (b) 
represents  the  fit  for  stations  0-17.  Correlation  coefficients  in  both 
cases  indicate  that  over  99%  of  the  variance  can  be  explained  by  a  linear 
fit. 
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Table  4.  CALIBRATION  COEFFICIENTS  FOR  REMOVAL  OF 
THE  SMALL  NEIL  BROWN  CTD  PRESSURE 
DEPENDENCE:  February  1989  has  two  separate  corrections 
based  upon  station  numbers.  November  1989  has  two  corrections 
based  upon  depth  range,  where  U  =  0-1000  dbar,  and  L  =  1000  dbar 
to  maximum  cast  depth. 


Cruise 

fa)  x-2 

(ai)x 

(ao) 

November  1988 

WKS& MIEMM 

-8.975E-6 

iKsom 

Heaassai i 

a— aaronai  n  hnnf 

MlillTM 

1.287E-5 

November  1989  -  U 

1.996E-5 

0.028 

November  1989  -  L 

0.0 

1.921E-6 

0.046 

Table  5.  SALINITY  CALIBRATION  COEFFICIENTS  FOR  THE 
SMALL  NEIL  BROWN  CTD:  The  standard  deviation  listed  in 
this  table  is  of  the  residuals  after  the  calibration  has  been  applied. 
The  number  of  points  indicates  the  number  of  instrument/bottle  pairs 
used  in  the  analysis. 


Cruise 

Slope 

Intercept 

#  of 
points 

Standard  Deviation 

November  1988 

1  -  17 

1.007040 

-0.246238 

28 

±0.012 

November  1988 
18-32 

1.002435 

-0.076711 

24 

±0.006 

February  1989 
0-17 

1.001894 

-0.065019 

21 

±0.011 

0.999418 

0.019968 

43 

±0.013 

November  1989 

1.002761 

-0.095870 

74 

±0.005 
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Once  the  pressure  effect  was  removed,  a  linear  fit  was  used  to  perform  the  final 
step  of  the  calibration.  The  coefficients  used  in  this  process  are  listed  in  Tables  4 
and  5. 

To  verify  the  effectiveness  of  these  calibrations,  a  deep  water  regime 
was  chosen  for  Temperature- Salinity  analysis.  Deep  water  masses  of  the  world 
oceans  have  been  regarded  as  having  relatively  stable  T-S  properties,  with  those 
of  the  North  Pacific  being  nearly  homogeneous  (Warren  and  Owen  1988).  Any 
variations  from  a  mean  T-S  curve  are  the  result  of  geographic  differences 
(Worthington  1981).  The  regime  selected  was  between  2800  and  3300  dbar. 
This  range  was  chosen  to  bracket  the  3000  dbar  level,  however,  the  limits 
themselves  were  arbitrary.  If  the  calibrations  were  successful,  the  deep  TS 
curves  for  cruises  along  the  POST  would  bear  a  strong  resemblance  to  one 
another,  and  predicted  salinity  values  would  approach  the  true  values  within  the 
accuracy  of  the  instrument.  After  plotting  these  deep  T-S  curves  (Appendix  A) 
for  each  cruise,  it  was  observed  that  in  fact  they  all  exhibited  a  similar  pattern 
and  were  all  within  ±  0.005  PSU  of  each  other.  From  these  results,  it  appeared 
the  calibration  procedures  were  successful  and  no  further  corrections  were 
necessary. 
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IV.  ANALYSIS 


A.  DESCRIPTION  OF  WATER  MASS  CHARACTERISTICS 

The  waters  off  Point  Sur  are  comprised  of  four  major  water  mass  types 
which  exist  in  varying  proportions,  summarized  in  Table  6.  The  location  and 
spatial  extent  of  the  water  masses  were  defined  using  an  analysis  similar  to  Lynn 
and  Simpson  (1990).  They  analyzed  the  flow  of  the  CUC  off  southern  California 
using  spiciness  and  spiciness  anomaly  as  tracers  to  identify  the  various  water 
masses  present  throughout  the  survey  region.  Spiciness  (71,  Flament  1986),  is  the 
state  variable  which  is  most  sensitive  to  isopycnal  thermocline  variations  and 
least  correlated  with  the  density  field.  Spiciness  is  useful  for  the  description  of 
interleaving  and  double  diffusive  processes  which  occur  at  the  boundary  between 
different  water  types.  Waters  which  are  warm  and  salty  have  positive  n  values 
while  those  which  are  cool  and  fresh  have  negative  rc  values. 

Table  6.  WATER  MASSES  OF  THE  CALIFORNIA  CURRENT 
SYSTEM:  (Source:  Simpson  1984) 


Water  Mass 

T 

s 

02 

Nutrients 

Surface  Water  Masses 
(0-200  m) 

Pacific  Subarctic 

Low 

Low 

■E391 

■iM 

North  Pacific  Central 

■n 

■EMM 

Low 

Coastally  Upwelled 

Low 

mmamm 

Low 

■m 

Subsurface  water  masses 

_ 

Equatorial  Pacific 

mnm* 

Low 

K333B 

Lynn  and  Simpson  (1990)  calculated  spiciness  anomalies  as  deviations  from  a 
mean  T-S  curve  representative  of  local  CC  water.  In  this  study,  only  profiles 
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free  of  interleaving  processes  below  the  surface  layer  (upper  50  m)  were  chosen 
for  the  average.  Eight  stations  from  the  seven  cruises  (Figure  11)  were  chosen 
for  this  average  (Figure  12)  and  were  assumed  representative  of  the  offshore 
waters  along  the  POST.  Anomalies  were  computed  along  level  surfaces  rather 
than  isopycnal  surfaces  for  ease  in  computation.  Anomalies  computed  along  level 
surfaces  only  differed  from  those  computed  along  isopycnal  surfaces  in  the  upper 
50  m  of  the  water  column.  Representative  curves  for  Pacific  Subarctic  and 
Equatorial  Pacific  waters  (Figure  12)  illustrate  that  mixing  has  produced  a 
profile  which  is  nearly  50%  Equatorial  Pacific  water  between  the  26.5  and  27.25 
density  surfaces  (roughly  250-750  m  depth)  and  predominantly  subarctic  above 
(Tibby  1941).  Because  the  calculation  of  spiciness  is  a  nonlinear  process,  the 
values  of  spiciness  for  each  of  the  eight  stations  were  computed  and  then 
averaged  (rather  than  vice-versa).  Due  to  the  exclusion  of  stations  where 
interleaving  processes  were  present  the  average  profile  in  the  upper  water 
column  was  dominated  by  those  stations  from  the  August  1988  cruise  where 
temperature  values  in  the  upper  water  column  were  high. 

1.  Analysis  of  spiciness  and  spiciness  anomaly 

Waters  of  subarctic  origin  (low  T  and  S)  will  be  represented  by 
negative  7t  values  in  the  Point  Sur  region  while  those  of  equatorial  origin  (high  T 
and  S)  will  be  represented  by  positive  n  values.  Due  to  the  strong  effects  of 
surface  heating,  the  study  of  spiciness  within  the  upper  50  m  was  difficult  and 
in  some  instances  had  to  be  neglected  in  subsequent  analysis.  The  following 
discussion  has  been  organized  by  cruises. 
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Figure  11.  The  eight  T-S  profiles  used  to  obtain  the  average  offshore 
profile.  Points  are  plotted  at  2  dbar  intervals. 


Figure  12.  Average  offshore  T-S  profile  used  in  calculation  of 
spiciness  anomalies:  Curves  representative  of  Pacific 
Subarctic  (long  dash)  and  Equatorial  Pacific  (short  dash)  waters 
from  Tibby  (1941)  are  superimposed  for  comparison. 
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a.  August  1988  Water  Mass  Analysis  (seasonal  normal) 

Data  from  the  August  1988  cruise  has  be  chosen  to  illustrate  the  relation  between 
the  fields  of  temperature,  salinity,  and  spiciness.  Inspection  of  the  T  and  S  fields 
(Figures  13  and  14,  respectively)  reveals  a  deepening  of  the  isotherms  and  a 
shallowing  of  the  isohalines  near  the  continental  slope.  This  is  characteristic  of 
the  CUC  which  carries  warmer,  saltier  water  poleward.  The  relatively  cooler 
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Figure  13.  Vertical  section  of  temperature  for  cruise  CUC-August  19S8: 

The  contour  interval  is  0.5°C. 
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Figure  14.  Vertical  section  of  salinity  for  cruise  CLC-August  1988: 

The  contour  interval  is  0.5  PSU. 

and  fresher  waters  of  Subarctic  origin  associated  with  the  CC  can  be  seen 
offshore  between  stations  16  and  22  below  200  m.  There  is  a  subsurface  salinity 
minimum  (50  m  depth),  with  values  less  than  32.8  PSU,  located  near  station  20 
associated  with  the  core  of  the  CC.  These  regions  are  reflected  in  the  vertical 
sections  of  spiciness  and  spiciness  anomaly  (Figures  15  and  16,  respectively). 
Comparison  of  Figures  15  and  16  reveals  a  positive  shift  in  magnitude  of 
roughly  0.2  units  in  spiciness  anomaly.  This  results  from  the  fact  that  anomalies 
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Figure  15.  Vertical  section  of  spiciness  for  cruise  CUC-August  1988: 

The  contour  interval  is  0.1  units. 


were  defined  as  deviations  from  the  average  offshore  profile.  !n  both  sections 
the  presence  of  the  CC  and  CUC  are  clearly  seen.  Negative  values  of  spiciness 
occur  between  stations  16  through  22  below  the  upper  50  to  75  m  with  a 
subsurface  minimum  (less  than  -0.2  units)  occurring  between  stations  20  and  21. 
This  corresponds  to  the  subsurface  salinity  minimum  discussed  earlier  and  is 
associated  witli  tlu  'ooler,  fresher  waters  of  the  CC.  At  station  20  there  is  also  a 
lensc  of  cooler  water  between  400  dbar  and  500  dbar  as  reflected  by  the  slopes 
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Figure  16.  Vertical  section  of  spiciness  anomaly  for  cruise  CUC-August 
1988:  The  contour  interval  is  0.1  units. 


of  the  5.5°C  and  6.0°C  isotherms  (Figure  13)  and  a  slight  depression  of  the 
isohalines  (Figure  14)  which  accounts  for  the  second  subsurface  minimum  (less 
than  -0.2  units).  It  is  believed  to  be  a  deeper  manifestation  of  the  Pacific 
Subarctic  waters  associated  with  the  CC.  Below  600  m,  negative  spiciness  values 
result  from  the  cooler  temperatures  naturally  occurring  with  depth.  Inshore 
(within  60  km  of  the  coast)  a  region  of  positive  spiciness  values  can  be  seen 
between  stations  1  and  15.  It  is  deepest  along  the  continental  slope  near  station  4. 
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This  region  of  positive  spiciness  is  associated  with  the  warmer,  more  saline 
waters  of  the  CUC,  and  lies  next  to  the  cooler  waters  of  subarctic  origin. 

Vertical  sections  of  spiciness  anomaly  offer  a  better  depiction  of  the 
spatial  extent  of  the  CUC,  and  further  discussion  will  be  limited  to  spiciness 
anomaly  although  mention  of  spiciness  will  be  made  when  necessary.  The 
spiciness  anomaly  section  for  August  1988  (Figure  16)  illustrates  the  same 
features  discussed  earlier,  however,  due  to  the  shift  in  magnitude  the  region  of 
positive  spiciness  anomaly  now  extends  to  a  depth  of  almost  700  m  along  the 
continental  slope.  This  agrees  favorably  with  the  position  of  the  CUC  seen  in  the 
August  1988  alongshore  geostrophic  velocity,  discussed  later  in  this  Chapter. 
Downward  sloping  isotherms  below  200  m  at  station  22  correspond  to  the  region 
of  positive  spiciness  anomaly  seen  in  Figure  16  and  are  suggestive  of  intrusion  of 
warm,  salty  water  which  entered  the  CC  from  the  west. 

b.  May  1988  Water  Mass  Analysis  ( weak  poleward  flow) 

A  period  of  strong  equatorward  wind  stress,  at  times  in  excess  of 
23  knots,  preceded  this  cruise  and  persisted  through  its  first  two  days.  These 
equatorward  winds  produced  upwelling  of  cold,  saline  water  near  the  coast  which 
was  then  transported  offshore  through  Ekman  dynamics.  Surface  fronts  were 
observed  in  the  vertical  fields  of  T  and  S  and  corresponded  to  those  seen  in  the 
density  field  (Figure  17). 

The  vertical  section  of  spiciness  anomaly  (Figure  18)  illustrates  a 
similar  pattern  with  positive  spiciness  anomaly  values  occurring  within  the  upper 
100-150  m  and  near  zero  anomalies  (of  Subarctic  origin)  below.  The  negative 
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Figure  17.  Vertical  section  of  density  anomaly  for  the  May  1988 
Student  Cruise:  The  contour  interval  is  0.1  units. 

values  found  near  the  surface  reflect  the  difference  from  the  average  profile 

which,  as  mentioned  earlier,  was  dominated  by  the  August  1988  cruise.  The 

region  of  positive  spiciness  anomaly  located  between  100  m  and  150  m  between 

stations  17  and  18  can  again  be  explained  as  an  intrusion  of  warm,  salty  water 

from  the  west.  There  is  a  region  of  positive  spiciness  anomaly  found  along  the 

continental  slope  centered  near  400  m  between  stations  2  and  8  which 

corresponds  to  an  area  of  weak  poleward  flow  associated  with  the  CUC.  The 
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relative  weakness  of  this  flow  may  be  related  to  the  strength  of  the  equatorward 
winds. 

c.  November  1988  Water  Mass  Analysis  ( strong  mesoscale 

* 

feature ) 

This  cruise  differed  from  the  rest  in  that  it  was  dominated  by  a 
rather  strong  anticyclonic  mesoscale  feature.  The  depression  in  the  isopycnals 
(Figure  19)  between  stations  4  and  14  have  the  characteristic  appearance  of  a 
deep  (greater  than  700  m)  warm-core  mesoscale  feature.  This  feature  is  also 
present  in  NOAA  9  AVHRR  satellite  imagery  (Figure  20)  taken  three  days  after 
the  occupation  of  POST  stations.  Breaker  and  Broenkow  (1989)  and  Tracy 
(1990)  discuss  a  recurring  anticyclonic  mesoscale  eddy  off  Monterey  Bay  whose 
position  varies  in  accordance  with  local  forcing.  It  is  hypothesized  that  the 
mesoscale  feature  observed  during  the  November  1988  cruise  is  a  southward 
displacement  of  this  recurring  anticyclonic  eddy.  Inshore  of  this  feature  the 
isopycnals  (greater  than  26.8)  slope  down  toward  the  continental  slope  in 
accordance  with  the  fields  of  temperature  and  salinity  (not  shown)  due  to  the 
presence  of  the  CUC. 

A  region  of  positive  spiciness  anomaly  (subsurface  maximum  near 
100  m)  is  associated  with  this  mesoscale  feature  indicative  of  relatively  warm, 
saline  water  characteristic  of  Equatorial  Pacific  origin  (Figure  21).  In  their 
study  of  a  mesoscale  dipole  eddy  off  southern  California,  Simpson  and  Lynn 
(1990)  found  that  anticyclonic  features  were  comprised  of  CUC  water.  Huyer  et 
al.  (1984)  analyzed  the  T-S  structure  of  a  mesoscale  feature  off  Oregon  and 
concluded  its  core  contained  water  of  shelf-slope  origin.  The  fact  that  larger 
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Figure  18.  Vertical  section  of  spiciness  anomaly  for  the  May  1988 
Student  Cruise:  The  contour  interval  is  0.1  units. 

values  of  spiciness  anomaly  (greater  than  0.3  unii  are  found  within  this  feature 

rather  than  along  the  continental  slope  suggests  that  lateral  entrainment  of  CUC 

water  has  occurred.  A  subsurface  minimum  (less  than  -0.5  units)  is  found 

between  stations  16  and  17  at  a  depth  of  50  m  corresponding  to  a  subsurface 

salinity  minimum  (less  than  32.9  PSU)  and  a  slight  doming  of  the  local  10.5  C 

isotherm  (not  shown).  This  is  characteristic  of  subarctic  waters  yet  occurs  in  a 

region  of  strong  poleward  (shown  later)  rather  than  equatorward  flow  as  would 

normally  be  expected. 
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Figure  20.  NOAA  AVHRR  satellite  image  from  0004  GMT, 

21  November  1988:  Darker  shades  along  the  coast  and  off 
Point  Sur  depict  colder  water.  Lighter  shades  found  off 
Monterey  Bay,  south  of  Point  Sur,  and  in  western  half  of  the 
image  depict  warmer  water.  Note  the  presence  of  an 
anticyclonic  feature  between  122°-  123°W,  and  36°-36.5°N. 
Grid  spacing  equals  30'. 
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Figure  21.  Vertical  section  of  spiciness  anomaly  for  cruise  CUC- 
November  1988:  The  contour  interval  is  0.1  units. 
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d.  February  1989  Water  Mass  Analysis  ( strong  equatorward 
flow  near  coast) 

Temperatures  within  the  upper  100  m  were  nearly  isothermal  with 
values  ranging  between  10.0°C  and  10.5°C.  This  agreed  favorably  with  NOAA 
1 1  AVHRR  satellite  imagery  (Figure  22)  for  the  cruise  period  which  also  depicts 


Figure  22.  NOAA  AVHRR  satellite  image  from  2058  GMT, 

5  February  1989:  Darker  shades  throughout  the  entire  image 
depict  colder  water  and  near  isothermal  conditions.  Grid  spacing 
equals  30'.  Note  the  nearly  isothermal  temperatures  between 
121.5°- 122.5°  W,  and  36°-36.5°  N. 
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an  isothermal  SST.  Within  45  km  of  the  coast,  isotherms  and  isohalines  (not 

shown)  below  100  m  exhibited  slight  downward  tilting  toward  the  continental 

slope.  In  the  surface  layer,  a  salinity  minimum  (less  than  33.3  PSU)  was  located 

0 

between  stations  5  and  6,  while  offshore  beyond  station  14  surface  salinity 
remained  between  33.5  PSU  and  33.6  PSU  throughout  the  upper  50  to  75  m. 

The  spiciness  anomaly  (Figure  23)  shows  values  greater  than  or 
equal  to  +0.1  over  the  range  from  75  m  to  500  m,  while  above  75  m  negative 
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values  prevail.  This  indicates  the  waters  below  75  m  were  relatively  warmer, 
and  saltier  that  the  average  offshore  profile,  and  may  result  from  the  strong 
CUC  which  reached  the  surface  (as  the  DC)  during  this  time. 

e.  May  1989  Water  Mass  Analysis  ( poleward  flow  near 

coast) 

A  vertical  section  of  salinity  (Figure  24)  depicts  several  dominant 
features.  First,  there  are  two  low  salinity  cores  (less  than  32.9  PSU)  between 
stations  12  and  13  and  at  station  17.  Second,  the  strong  halocline  centered  at  150 
m  offshore  shoals  to  less  than  50  m  inshore.  Last,  there  are  large  amplitude 
variations  in  the  34.1  and  34.2  isohalines.  Along  the  continental  slope,  isotherms 
slope  downward  toward  shore  (not  shown)  and  isohalines  upward  (Figure  24) 
toward  shore  indicating  the  water  is  warmer,  and  more  saline  than  offshore. 
The  depressions  in  the  the  34.1  isohaline  and  the  low  salinity  cores  correspond  to 
regions  of  equatorward  flow  (shown  later)  and  likely  contain  Pacific  Subarctic 
water  associated  with  the  CC. 

The  cores  of  low  salinity  cause  the  regions  of  negative  spiciness 
anomaly  above  100  m  depth  (Figure  25).  Below  200  m,  waters  with  negative 
spiciness  values  of  Subarctic  origin  appear  as  a  blank  region  below  an  overlying 
region  of  positive  anomaly  values,  comprised  of  two  different  water  masses. 
The  first  occurs  along  the  continental  slope  where  values  greater  than  +0.2  units 
are  observed  and  values  of  +0.1  extend  downward  to  600  m.  This  region  is 
associated  with  the  Equatorial  waters  of  the  CUC.  Offshore  the  +0.1  contour  has 
a  subsurface  tongue-like  appearance  extending  to  approximately  240  km  from 
shore  at  150  m  depth  near  stations  15  and  16.  A  subsurface  maximum 
(greater  than  +0.2  units)  is  found  between  stations  1 1  and  12  corresponding  to  a 
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subsurface  anticyclonic  feature  present  in  the  density  field  (Figure  26).  As  in  the 
case  of  the  mesoscale  feature  observed  in  November  1988,  it  is  hypothesized 
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Figure  24.  Vertical  section  of  salinity  for  the  May  1989  Student  Cruise: 

The  contour  interval  is  0.1  PSU. 


that  these  larger  values  offshore  result  from  lateral  entrainment  and  offshore 
advcction  of  CUC  w  ater. 

Further  offshore  near  station  19  is  another  region  of  positive 
spiciness  anomaly  believed  to  result  from  the  intrusion  of  warmer,  saline  water 
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from  the  west.  This  intrusion  can  be  seen  along  the  24.75  isopycnal  in  the  T-S 
curve  for  station  19  (Figure  11)  and  will  be  discussed  further  in  Chapter  V. 


Figure  25.  Vertical  section  of  spiciness  anomaly  for  the  May  1989 
Student  Cruise:  The  contour  interval  is  0.1  units 

/.  July  1989  Water  Mass  Analysis  (anomalous  poleward 
flow) 

The  water  mass  structure  shown  in  the  vertical  section  of  spiciness 
anomaly  (Figure  27)  is  very  different  than  August  1988  (Figure  16).  one  year 
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earlier.  The  section  is  almost  completely  dominated  by  positive  values  of 
spiciness  anomaly,  which  extend  to  depths  of  600  m  in  some  locations.  Prior  to 
and  during  the  cruise  the  winds  at  both  buoy  locations  were  equatorward 


Figure  26.  Vertical  section  of  density  anomaly  for  the  May  1989 
Student  Cruise:  The  contour  interval  is  0.1  kg  m  V 

(<  -1.0  dyne  enr2),  howeve  a  substantial  wind  relaxation  lasting  3  days  preceded 

the  cruise  by  two  days.  Remote  forcing  by  this  wind  relaxation  is  believed  to  be 

one  of  the  mechanisms  involved  in  generating  this  anomalous  field.  Surface 

salinities  along  the  section  on  this  cruise  (Figure  28)  are  approximately  0.6  PSU 
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greater  than  the  August  1988  cruise.  Further  discussion  of  this  relaxation  and 
the  associated  poleward  flow  will  be  presented  later  in  this  Chapter. 


Figure  27.  Vertical  section  of  spiciness  anomaly  for  cruise  CUC-July 
1989:  The  contour  interval  is  0.1  units. 


T-S  curves  for  stations  20  and  21  (not  shown)  reflect  the  intrusion 
of  warmer,  saline  water  along  the  25.1  kg  nv3  isopycnal  surface.  Water  inshore 
of  stations  18  and  19  are  composed  of  Equatorial  Pacific  water  advected 
northward  with  the  observed  strong  poleward  flow'.  The  core  of  positive 
anomaly  (>  +0.3  units)  is  associated  with  the  CUC  which  appears  to  have 
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surfaced  during  this  cruise.  Below  150  m  there  is  a  region  of  negative  spiciness 
anomaly  occurring  between  stations  18  and  22  believed  to  be  of 


Figure  28.  Vertical  section  of  salinity  for  cruise  CUC-Julv  1989: 

The  contour  interval  is  0.1  PSU. 

Subarctic  origin  associated  with  the  CC,  even  though  it  occurs  in  a  region  of 
poleward  flow  (>  5  cm  s*1,  shown  later).  Centered  near  a  depth  of  50  m  between 
stations  17  and  19  are  small  areas  of  negative  anomalies  associated  with  CC 
water.  Inshore  between  stations  11  and  13  is  a  core  of  negative  anomaly  values 
associated  with  a  core  of  low  salinity  water  (Figure  28).  T-S  curves  for  stations 
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11  through  13  (not  shown)  reflect  the  intrusion  of  cool,  low  salinity  water  along 
the  26.8  kg  nr3  isopycnal  surface  (depth  approximately  50  m),  characteristic  of 
Subarctic  origin.  This  core  lies  within  a  region  of  poleward  flow  and  is 
probably  the  northward  transport  of  CC  water  associated  with  the  observed  wind 
relaxation. 

g.  November  1989  Water  Mass  Analysis  ( deep  equatorward 

flow) 

November  1989  was  the  final  cruise  included  in  this  study  and 
shows  marked  differences  from  the  November  1988  cruise,  one  year  earlier. 
Isopycnal  surfaces  below  250  m  (Figure  29)  have  a  domed  appearance  which 
extends  over  the  entire  transect.  They  are  closest  to  the  surface  near  station  19 
and  slope  downward  away  from  this  station.  Above  250  m,  isopycnals  between 
stations  21  and  22  begin  to  slope  upward  as  they  progress  offshore,  while  those 
inshore  of  station  19  continue  to  slope  up  toward  the  coast  until  they  reach 
stations  7  and  8  where  they  begin  to  slope  down  toward  the  continental  slope.  A 
sharp  density  front  is  located  near  station  17  and  marks  a  transition  region 
between  equatorward  flow  inshore  and  poleward  flow  offshore  (shown  later). 
The  corresponding  vertical  section  of  spiciness  anomaly  (Figure  30)  shows 
several  interesting  features.  First,  the  upper  100  m  of  the  transect  is 
characterized  by  a  complex  series  of  sign  reversals  in  the  anomaly  field.  Areas 
of  subarctic  waters  are  shown  by  the  negative  values,  while  those  of  southern 
origin  are  shown  by  positive  values.  Inshore  the  deeper  signature  of  the  CUC  is 
seen,  with  a  subsurface  maximum  (>  +0.3  units)  lying  between  stations  3  and  7  at 
a  depth  of  approximately  110  m.  Positive  values  extend  downward  to  a  depth  of 
500  m,  which  agrees  favorably  with  the  poleward  flow  of  the  CUC  (shown 


56 


later).  Second,  a  region  of  negative  spiciness  values  within  the  interior  of  the 
transect,  shown  as  the  large  blank  area  centered  between  stations  12  and  21 


CUC-NOVEMBER  1989  1 — 1 

DENSITY  ANOM  (?)  10  KM 


Figure  29.  Vertical  section  of  density  anomaly  for  cruise  CUC-November 
19K9:  The  contour  interval  is  0.1  kg  m-T 

(Figure  30).  This  region  corresponds  to  the  doming  of  the  isotherms  (not 
shown)  and  depression  of  the  isohalines  (Figure  31),  which  are  reflected  in 
the  density  field  (Figure  29).  Lastly,  the  extensive  region  of  positive  anomalies 
found  between  stations  21  and  22  below  100  m  depth  results  from  the  downward 
sloping  isotherms  and  sharply  upward  sloping  isohalines  resulting  in  water  which 
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is  warmer,  and  more  saline  than  adjacent  inshore  waters.  This  water  may  either 
be  a  very  strong  and  deep  manifestation  of  North  Pacific  Central  water  entering 
the  CCS  from  the  west  at  depth,  or  an  unusually  strong  recirculation  of  CUC 
water.  North  Pacific  Central  water  has  a  salinity  minimum  of  approximately 
34.1  PSU  between  depths  of  300  m  to  800  m  (Pickard  and  Emery  1982).  The 
resulting  slope  of  the  isopycnals  has  produced  extremely  deep  equatorward  flow 
(flows  <  -5.0  cm  s-1  at  700  m). 


November  1989:  The  contour  interval  is  0.1  units. 
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Figure  31.  Vertical  section  of  salinity  for  cruise  CUC-November  1989: 

The  contour  interval  is  0.1  PSU. 


B.  DESCRIPTION  OF  ALONGSHORE  GEOSTROPHIC  FLOW 
1.  Geostrophy  and  its  limitations 

In  the  absence  of  wind  stress  and  frictional  forces,  along  with  proper 
scaling  considerations,  the  equations  of  motion  (la,b,c)  can  be  reduced  to  (2a, b) 
in  the  horizontal  and  (2c)  in  the  vertical.  The  former  are  commonly  referred  to 
as  the  geostrophic  equations  while  the  latter  is  known  as  the  hydrostatic  equation. 
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These  equations  are  suitable  for  large  scale  ocean  currents,  where  the 
flow  is  assumed  to  be  steady  and  the  primary  balance  is  between  the  pressure 
gradient  and  Coriolis  forces.  This  involves  phenomena  with  length  scales  longer 
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than  the  first  internal  Rossby  radius  of  deformation.  To  compute  relative 
geostrophic  velocities  from  an  observed  density  field  we  need  to  redefine  these 
equations  in  terms  of  density  instead  of  pressure.  This  ca/i  be  accomplished  by 
taking  the  vertical  derivatives  of  equations  (2a,b)  and  the  horizontal  derivative  of 
equation  (2c).  Combing  these  results  we  obtain  what  are  commonly  referred  to 
as  the  "thermal  wind"  equations: 


3v  _  g  9p 

dz  p  fdx  (3a) 

du  _  g  dp 

dz  pfdy  (3b) 

From  these  two  equations  it  is  possible  to  compute  the  vertical  shear  of  the 
geostrophic  velocity  provided  the  horizontal  gradient  of  density  is  known.  In 
practice,  however,  it  is  desirable  to  use  a  variation  of  these  equations,  which  is 
based  upon  the  geopotential  distance  between  two  pressure  surfaces. 

The  "geostrophic  method"  as  described  by  Pond  and  Pickard  (1983) 
utilizes  the  temperature  and  salinity  information  collected  at  each  station  to 
compute  the  specific  volume  anomaly,  8.  The  specific  volume  anomaly  is  the 
sample  specific  volume,  a,  minus  the  specific  volume  of  standard  sea  water  (8  = 
c*s,t,p  -  cx35  o  p).  Using  this  quantity,  the  geostrophic  shear  (m  s'1)  between  two 
pressure  surfaces  can  be  computed  from  data  at  stations  A  and  B  as  (Pond  and 
Pickard  1983): 
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v2  -  Vi  = - i - [ADb  -  AD  a] 

L2Qsin<t> 


is  the  dynamic  height  anomaly  (m2  s-2)  at  each  station,  L  the  station  spacing 
(meters),  5  in  m3  kg1,  and  p  in  pascals. 

The  velocity  resulting  from  this  equation  is  the  relative  shear  between 
two  arbitrary  isobaric  surfaces.  To  convert  this  relative  velocity  to  an  absolute 
velocity  we  may  consider  two  possibilities.  First,  a  inference  "level  of  no 
motion",  vref  =  0,  can  be  assumed  where  the  corresponding  isopycnal  surface 
between  station  pairs  is  assumed  to  be  level.  Second,  additional  information 
from  a  current  meter,  ADCP,  or  PEGASUS  instrument  can  be  used  to  prescribe 
a  "level  of  known  motion",  where  vref  is  not  necessarily  equal  to  zero.  The 
geostrophic  velocities  discussed  throughout  the  remainder  of  this  paper  have 
been  computed  using  the  assumption  of  a  level  of  no  motion  at  1000  dbar  based 
upon  the  PEGASUS  observations  of  Rago  and  Collins  (1989)  along  the  POST 
from  April  1988  to  March  1989. 

Before  proceeding  further  it  is  appropriate  to  first  discuss  the 
limitations  of  (or  errors  in)  the  geostrophic  method  which  arise  due  to  the 
assumptions  that  have  been  made  in  its  calculation.  The  geostroplrc  method  has 
several  disadvantages  (Pond  and  Pickard  1983)  which  are  as  follows; 

-  it  produces  a  relative  current  which  itself  is  dependent  upon  the  assumption 
of  a  level  of  no  motion; 
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-  a  problem  arises  when  the  selected  reference  level  reaches  or  becomes 
deeper  than  the  ocean  bottom  which  occurs  as  stations  become  close  the 
shore; 

-  it  is  best  suited  for  stations  which  are  tens  of  kilometers  apart  where  it  will 
yield  a  mean  velocity  over  that  distance; 

-  friction,  which  has  been  neglected,  may  in  fact  be  large  near  the  bottom  or  in 
regions  where  there  is  strong  shear; 

-  the  equations  do  not  apply  near  the  equator  (±0.5°)  where  the  Coriolis 
parameter  tends  toward  zero; 

-  internal  waves  and  tides  can  disturb  the  mass  field  from  its  equilibrium  state 
aliasing  the  data  in  an  unknown  way. 

2.  Calculation  of  alongshore  geostrophic  velocity 
Geostrophic  velocities  have  been  computed  based  upon  a  reference  level 

of  1000  dbar.  The  method  used  here  to  extrapolate  nearshore  is  that  utilized  by 
Reid  and  Mantyla  (1976),  in  which  the  value  of  dynamic  heights  have  been 
linearly  extrapolated  horizontally  along  their  last  observed  slope  (Figure  32). 
This  technique,  though  artificial,  has  been  used  by  many  investigators  (Huyer 
1980;  Lynn  and  Simpson  1987;  Lynn  and  Simpson  1989)  in  the  calculation  of 
alongshore  geostrophic  velocities  and  appears  to  be  the  preferred  technique. 

3.  Estimation  of  error  in  dynamic  height  and  geostrophic 
velocity 

Besides  the  inherent  limitations  of  geostrophy,  the  observations 
themselves  often  contain  error,  which  has  a  direct  effect  on  the  final  velocity 
field.  Throughout  this  paper  the  term  "error"  refers  to  random  error,  and 
"bias"  to  systematic  errors.  Random  error  is  equivalent  to  the  precision  of  the 
quantity  in  question  and  is  not  a  measure  of  accuracy.  Accuracy  includes  both 
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Figure  32.  Extrapolation  of  dynamic  height  relative  to  1000 
decibars:  This  Figure  illustrates  the  technique  of  a  linear 
extrapolation  of  the  dynamic  height  field  to  estimate  geostrophic 
currents  when  the  bottom  becomes  shoaler  than  the  reference 
level.  The  technique  is  based  upon  using  the  last  observed  slope 
of  dynamic  height  to  project  what  the  inshore  values  might  be  if 
the  bottom  were  replaced  by  a  Fictitious  water  column.  The 
dashed  line  at  the  surface  beginning  at  station  65  represents  the  sea 
surface  value  which  would  result  if  the  dynamic  heights  were 
extended  horizontally  instead  of  along  the  last  observed  slope. 
(Source:  Reid  and  Mantyla  1976) 
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random  errors  and  bias.  Random  error  is  expressed  as  the  standard  deviation  of 
the  errors  (o),  and  a2  is  the  variance  of  the  distribution. 

Reid  (1958)  discusses  several  sources  of  error  .which  can  effect  the 
dynamic  computations.  When  necessary  this  discussion  has  concentrated  on  the 
more  modern  instrumentation  used  today.  Possible  sources  include; 

-  incorrect  temperatures  which  arise  from  improper  calibration  of  the  CTD 
temperature  sensor; 

-  incorrect  salinities  which  arise  from  either  contami  ion  of  water  samples, 
operator  error  on  the  lab  salinometer,  or  from  errors  in  the  calibration 
procedure; 

-  incorrect  estimate  of  the  pressure  at  which  the  readings  and  samples  were 
taken  due  to  improper  calibration  of  the  CTD  pressure  sensor; 

-  improper  station  location  and  spacing  due  to  errors  in  navigation  and  vessel 
drift; 

-  errors  in  the  mass  field  due  to  inadequate  sampling  which  can  result  from 
irregularities  in  the  horizontal  and  vertical  mass  structure,  fluctuations  of  the 
field  with  time,  and  the  fact  that  stations  are  not  occupied  simultaneously; 

-  and  displacement  of  isopycnals  due  to  internal  waves. 


Both  Reid  (1958)  and  Fofonoff  (1985)  discuss  the  effect  that  errors  in 
the  coefficients  of  the  equation  of  state  would  have  on  the  dynamical 
computations  described  earlier.  Although  different  versions  of  this  equation 
were  used,  both  authors  conclude  that  such  errors  are  insignificant  when 
compared  to  those  which  result  from  internal  tides  and  measurement  errors.  For 
this  reason  the  errors  in  the  coefficients  of  the  equation  of  state  will  be  neglected 
in  subsequent  error  analysis  of  geostrophic  velocities. 
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To  determine  the  effect  of  errors  in  the  measurement  of  temperature, 
salinity  and  pressure  on  the  computation  of  dynamic  height  and  geostrophic 
velocity  an  analysis  similar  to  that  of  Wooster  and  Taft  (1958)  and  Johns  (1984) 
has  been  performed.  Briefly,  this  process  consists  of:  1)  propagating  the  errors 
of  T,  S,  and  P  through  the  equation  of  state  to  determine  the  corresponding  error 
in  the  specific  volume  anomaly;  2)  translating  this  into  an  error  in  dynamic 
height;  and  3)  using  this  value  along  with  the  error  in  station  spacing  to 
determine  the  error  in  geostrophic  velocity.  Each  step  in  this  procedure  is 
described  next  in  detail. 

a .  Determination  of  error  in  specific  volume  anomaly 

To  propagate  the  measurement  errors  into  the  specific  volume 
anomaly  it  is  necessary  to  differentiate  the  equation  of  state  with  respect  to 
temperature,  salinity  and  pressure.  The  equation  of  state  used  here  is  that  found 
in  UNESCO  Technical  Report  #44  and  Millero  and  Poisson  (1981).  Since 
analytic  differentiation  of  this  equation  proves  to  be  a  difficult  task,  these 
derivatives  have  been  computed  numerically  through  finite  differences. 

Twenty  values  of  temperature,  salinity  and  pressure  were  randomly 
chosen  from  four  separate  cruises  to  use  the  maximum  and  minimum  values  of 
these  variables.  Keeping  two  of  the  three  variables  fixed  the  third  was  increased 
by  either  0.005°C,  0.005  PSU,  or  1%  of  depth.  These  values  represent  the 
accuracies  of  the  Neil  Brown  CTD,  and  are  subsequently  used  as  estimates  of  the 
la  random  error  in  the  measurements.  The  resulting  values  for  these  derivatives 
were; 


jS  =  0.2  x  IQ'6 
dT  °C 
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—  =  -0.7  x  10‘6 

as  psu 

r 

—  =  0.5  x  10~'om3^-' 

9P  dbar 

It  is  clear  that  the  dominant  terms  are  those  associated  with  temperature  and 
salinity  which  are  four  orders  of  magnitude  largei  than  the  pressure  term. 

Using  these  values  it  is  possibk  compute  the  corresponding  error 
in  the  specific  volume  anomaly  in  a  manner  similar  to  thaf  of  Wooster  and  Taft 
(1958).  If  it  is  assumed  that  the  measurements  are  statistically  independent  the 
following  expression  for  the  variance  of  specific  volume  anomaly  (Gg2)  can  be 
utilized; 
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Substitution  of  these  derivatives  along  with  the  estimated  precision  of  the  Neil 
Brown  CTD  into  this  equation  yields  a  variance  of  1.325  x  10  ",7  I'm7  kg1)2 
which  corresponds  to  an  ‘’rror  of  ±  3.64  x  10  '9  m3  kg-1.  The  value  used  for  op2 
was  2  dbar.  Development  of  this  equation  is  based  upon  the  principles  outline;’ 
in  Appendix  t>.  The  equation  used  by  Wooster  and  Taft  (1958)  in  their  analysis 
included  an  additional  term; 


oj-  =  Oc2  +  ( 


(6) 
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where  the  first  term  represents  the  variance  of  the  measurement  errors, 
analogous  to  the  equation  shown  earlier,  and  the  second  term  represents  the 
variance  of  depth  errors.  They  find  that  below  the  surface  layer  this  second 
term  is  much  smaller  than  the  first,  even  if  large  vertical  gradients  are  present  at 
depth,  and  can  therefore  be  neglected  in  the  computation  of  dynamic  height 
errors.  Further,  they  indicate  that  the  first  term  is  dominated  by  salinity,  which 
agrees  with  the  results  of  this  analysis.  For  the  bulk  of  this  analysis  this  term  has 
been  neglected,  however  it  will  be  included  when  the  effects  of  internal  tides  are 
examined. 

b.  Determination  of  error  in  dynamic  height 

The  next  step  is  to  translate  the  error  in  the  specific  volume 
anomaly  into  an  error  in  dynamic  height.  In  the  process  of  determining  the 
derivative  of  5  with  respect  to  temperature,  salinity  and  pressure,  it  was  found 
that  the  measurement  errors  did  not  vary  significantly  with  these  variables.  As  a 
result,  the  calculation  of  the  error  in  dynamic  height  will  depend  on  the  number 
of  levels  used  in  the  integration.  In  calculating  the  geostrophic  velocities,  a  layer 
thickness  of  20  dbar  has  been  used.  Since  the  data  has  been  averaged  into  2  dbar 
bins,  the  error  in  dynamic  height  must  first  be  computed  within  this  layer,  and 
then  carried  from  the  reference  level  to  the  surface  and  bottom.  The  variance  in 
dynamic  height  within  a  20  dbar  layer  has  been  computed  from  the  following 
equation; 
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where  pressure  has  units  of  Pascals  (kg  m  s'2  m*2).  Using  the  value  obtained 
earlier  for  the  error  in  specific  volume  anomaly,  the  resulting  error  in  dynamic 
height  for  a  20  dbar  layer  becomes  ±  0.00023  m  2  s  '2.  Since  the  error  is  the 
same  at  each  station  in  the  pair  the  value  used  to  calculate  the  error  in  the 
geostrophic  velocity  is  simply  the  square-root  of  the  sum  of  the  squares  which  is 
±  0.00033  m2  s'2. 

To  obtain  the  random  error  at  any  surface  relative  to  1000  dbar, 
multiply  this  value  by  the  square-root  of  the  number  of  layers.  For  the  sea 
surface  relative  to  1000  dbar  (50  layers)  a  random  error  in  dynamic  height  of  ± 
0.0023  m2  S'2  or  ±  0.00023  dynamic  meters  is  estimated  as  error  in  height 
difference  between  two  stations.  The  smallness  of  this  result  can  be  attributed  to 
good  measurements  and  to  the  fact  that  the  data  was  continuous  which  allowed 
for  the  use  of  small  layers  in  the  computational  process.  This  in  turn  allowed  for 
a  better  estimate  of  the  integral  in  equation  (4). 

This  value  is  similar  to  that  obtained  by  Johns  (1984),  ±  0.0006 
dynamic  meters  relative  to  2000  dbar,  where  the  small  size  results  from  the  fact 
that  the  measure  its  were  assumed  to  be  free  of  any  biases.  As  she  points  out, 
this  value  is  very'  small  and  probably  underestimates  the  total  (random  plus  bias) 
error  in  the  dynamic  height.  To  obtain  a  more  realistic  result  she  allows  for  a 
bias  of  ±  0.0025  °C  and  ±  0.0025  PSU  in  the  measurements  of  temperature  and 
salinity  and  obtains  a  value  of  ±  0.004  dynamic  meters. 

Application  of  biases  in  the  temperature  and  salinity  measurements 
will  certainly  have  an  effect  on  the  value  of  dynamic  height  and  will  be 
important  when  considering  plots  of  dynamic  topography.  However,  in  the 
computation  of  geostrophic  velocities  it  is  the  difference  in  dynamic  heights 


between  two  stations  which  is  important  and  a  fixed  bias  will  tend  to  cancel  out. 
For  this  reason  no  bias  has  been  introduced  into  this  analysis. 

c.  Determination  of  error  in  geostrophic  .velocity 

To  obtain  an  equation  for  the  error  in  geostrophic  velocity  take  the 
total  derivative  of; 


vg  =  (ADb  -  ADa)  +  vrcf 

IL  (6) 

and  the  square-root  of  the  sum  the  squares  of  the  components  to  obtain  (assuming 
errors  in  L,  vrCf,  and  AD  are  statistically  independent); 


where  L  is  in  meters,  vg  is  in  m  s  *J,  Oad  is  in  m2  S'2,  and  Ol  is  in  meters 
(Appendix  B).  These  three  terms  represent  the  error  in  the  dynamic  height,  the 
error  in  station  spacing  and  the  error  in  the  velocity  at  the  reference  level, 
respectively.  Prior  to  computing  error  estimates  for  geostrophic  velocities,  the 
error  in  station  spacing  must  First  be  determined. 

Unlike  previous  error  quantities  which  do  not  vary  significantly 
from  cruise  to  cruise,  the  error  in  station  spacing  must  be  calculated  separately 
for  each  cruise.  The  two  sources  of  error  in  station  spacing  are  navigational 
errors  and  vessel  drift.  The  navigation  error  is  insignificant  compared  to  the 
vessel  drift  which  depends  upon  weather  conditions  and  speed  of  surface 
currents.  The  navigation  system  used  to  position  the  vessel  was  LORAN-C, 
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which  for  the  coastal  waters  in  the  vicinity  of  Point  Sur  has  an  absolute  accuracy 
of  about  0.25  nautical  miles  (Bowditch  1984).  Most  of  this  error  is  a  spatially 
coherent  bias  which  cancels  upon  differencing  station  positions  to  determine  the 
station  spacing.  The  remaining  random  error  in  station  spacing  is  on  the  order 
of  approximately  ±  25  m  (K.J.  Schnebele,  pers.  comm.)  in  this  region  of  the 
coast.  This  is  a  very  small  and  essentially  negligible  in  the  geostrophic 
calculations.  Therefore  it  is  necessary  to  compute  the  differential  drift  between 
all  applicable  stations  pairs  for  each  cruise  to  determine  an  appropriate  value  for 
Ol  to  be  used  in  the  geostrophic  velocity  error  calculations. 

From  the  starting  and  ending  positions  for  each  cast,  the  component 
of  the  change  in  station  spacing  which  lies  along  the  POST  has  been  calculated 
for  each  pair.  If  at  both  stations  the  vessel  drifted  in  the  same  direction  the 
values  were  subtracted,  while  if  the  drift  was  in  opposite  directions  the  values 
were  added.  Since  continuous  data  was  collected  on  the  downcast  the  maximum 
uncertainty  in  the  spacing  is  one  half  this  distance.  This  assumes  that  the  ship 
drifted  one-half  the  total  distance  during  the  data  collection  process  and  is  valid 
for  all  casts,  except  those  where  numerous  bottle  samples  were  collected  on  the 
upcast.  In  these  instances  the  actual  drift  would  be  less  than  one-half  so  that  this 
procedure  would  tend  to  produce  more  conservative  results.  If  it  is  also  assumed 
that  the  drift  of  the  instrument  on  the  downcast  was  linear  with  time  the  error  is 
zero  at  the  surface  and  reaches  the  maximum  uncertainty  at  the  bottom. 

Based  upon  the  uniform  probability  distribution  an  rms  value  was 
determined  by  dividing  the  maximum  uncertainty  by  2  times  the  square-root  of  3 
(Appendix  B).  This  value  was  then  divided  by  the  station  spacing  to  obtain  the 
dimensionless  parameter,  Ol  /  L,  needed  for  the  calculation  of  the  error  in 
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velocity.  For  each  cruise  an  average  value  has  been  obtained,  along  with  the 
maximum  and  minimum  values  (Table  7). 

The  third  term  in  this  equation  represents  tjie  uncertainty  of  the 
flow  at  the  reference  level.  Based  upon  PEGASUS  observations  of  Rago  and 
Collins  (1989),  the  flow  at  1000  dbar  is  near  zero  ±  0.02  ms*1.  Where  the 
bottom  depth  at  a  station  was  less  than  the  reference  level  extrapolation  of 
dynamic  height  was  required.  Since  this  extrapolation  was  linear,  the  deepest 
velocity  in  the  extrapolated  pair  was  equal  to  the  velocity  of  the  adjacent  offshore 
pair  at  the  same  depth.  Therefore  the  random  error  at  the  deepest  common 
depth  of  the  next  inshore  station  pair  is  equal  to  the  value  of  avg  of  the  adjacent 
offshore  station  pair. 

Table  7.  NORMALIZED  ALONGTRACK  DIFFERENTIAL  DRIFT 
OF  THE  VESSEL:  These  values  represent  the  rms  differential 

drift  normalized  by  the  spacing  between  station  pairs  (  ol  /  L ).  Drift 
refers  to  a  distance  rather  than  a  measure  of  current  speed. 


Cruise 

Average 

Maximum 

Minimum 

STMAY  1988 

0.0088 

0.025 

0.0001 

mm 

0.0157 

0.059 

0.0004 

November  1988 

0.0207 

0.062 

0.001 

February  1989 

0.0155 

0.055 

0.000 

STMAY  1 989 

0.0114 

0.022 

0.001 

■Hi 

0.0154 

0.041 

0.0003 

November  1989 

0.0170 

0.044 

0.002 
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Using  the  maximum  value  for  Ol  /  L  from  (Table  7),  the  value  for  Gad 
determined  earlier,  and  the  appropriate  value  for  the  Coriolis  parameter,  /, 
which  for  this  latitude  is  8.6408  x  10  5,  the  worst  case  error  in  the  geostrophic 
velocity  can  now  be  determined.  The  November  1988  cruise  had  the  highest 
values  for  both  ctl  /  L  (0.062)  and  geostrophic  velocities  (45  to  55  cm  s*1)  and 
was  used  for  a  test  case  to  determine  the  magnitude  of  the  resulting  errors  in 
geostrophic  velocities.  The  errors  were  ±  2.0  cm  s'1  at  the  reference  level  and 
increased  to  a  maximum  near  the  surface  or  where  the  geostrophic  currents  were 
largest  and  at  the  deepest  depth  below  the  reference  level.  For  this  particular 
cruise  the  largest  errors  were  on  the  order  of  approximately  3.5  to  4.0  cm  s_1, 
which  occurred  in  the  areas  of  the  greatest  geostrophic  velocities.  Based  upon 
these  results  it  appears  that  the  worst  case  error  in  station  spacing  contributes  to 
an  error  of  approximately  ±  1.5  -  2.0  cm  s'1.  The  errors  for  other  cruises  were 
much  smaller,  consisting  primarily  of  the  ±  2.0  cm  s1  error  in  the  reference 
level  velocity. 

d.  Effect  of  internal  tides  on  errors  in  dynamic  height 

Baroclinic  tides  result  from  the  interaction  of  barotropic  tides  with 
variable  bottom  topography  associated  with  coastal  regions.  Vertical  particle 
displacements  of  10  to  100  meters  and  horizontal  velocities  of  0.10  to  1.0  m/sec 
can  result  from  internal  tides  (Wiseman  et.al.  1983).  Observations  of  internal 
tides  over  six  days  off  Point  Sur,  California  in  October  1950  (Reid  1956),  show 
wave  amplitudes  greater  than  30  feet  (9.1  m). 

To  examine  the  effect  of  internal  tides  on  dynamic  height  and 
geostrophic  velocities  it  is  necessary  to  return  to  equation  (6)  (Wooster  and  Taft 
1958)  and  include  the  second  term.  The  first  part  of  this  term  represents  a 
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typical  vertical  gradient  of  the  specific  volume  anomaly,  8,  while  the  second  part 
represents  the  variance  associated  with  the  sinusoidal  internal  wave.  A  typical 
value  for  the  gradient  of  specific  volume  anomaly,  6.92381  x  10  9  kg  nr3  /  nr1, 
was  obtained  by  observing  the  change  of  8  over  a  20  meter  distance  in  the 
vicinity  of  the  thermocline  for  the  May  1988  student  cruise.  For  an  internal  tide 
wave  with  amplitude  10  m  the  variance  becomes  50  m2.  Combining  these  values 
with  the  value  of  the  first  term  and  following  the  same  procedure,  a  value  of  ± 
0.0031  m2  S'2  is  obtained  for  the  random  error  in  dynamic  height  of  a  20  dbar 
layer  at  each  station.  Notice  that  this  value  is  an  order  of  magnitude  larger  than 
that  obtained  without  the  inclusion  of  an  internal  tide.  It  should  be  kept  in  mind 
that  this  value  occurs  in  the  thermocline  region  where  the  value  of  the  specific 
volume  anomaly  changes  most  rapidly.  Above  and  below  this  region  the  rate  of 
change  and  thus  the  effect  of  the  internal  tide  is  much  smaller. 

Defant  (1950)  discusses  the  difficulties  associated  with  internal  tides 
and  the  resulting  periodic  fluctuations  in  the  mass  field.  From  observations  off 
southern  California,  he  found  that  internal  tides  can  cause  disturbances  of  up  to  8 
dynamic  centimeters.  He  also  discusses  a  method  for  removing  the  effect  of  the 
diurnal  or  semidiurnal  tides  through  careful  preparation  of  the  sampling  scheme. 
This  procedure  although  capable  of  removing  a  desired  tidal  frequency  is  not 
very  practical  from  an  operational  point  of  view.  In  addition,  the  period  of  the 
local  internal  tide  is  not  known  in  most  cases  and  is  therefore  extremely  difficult 
if  not  impossible  to  account  for  or  remove  its  effects  on  the  mass  field.  As  a 
result  the  reader  should  be  aware  that  internal  tides  do  exist  and  will  have  an 
effect  on  the  observed  mass  field  and  resulting  dynamic  topography,  however. 
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without  actual  measurements  of  these  tides  their  effects  cannot  be  isolated  and/or 
removed. 

4.  Analysis  of  alongshore  geostrophic  velocity 

* 

A  description  of  alongshore  geostrophic  velocities  for  each  of  the  seven 
cruises  is  presented.  With  the  exception  of  the  November  1989  California 
undercurrent  cruise  a  reference  level  of  1000  dbar  appeared  justified. 
Comparisons  with  earlier  studies  off  Point  Sur  and  discussion  of  possible 
dynamics  involved  follow  in  Chapter  V.  Observations  are  presented  by  cruises 
following  a  format  similar  to  the  water  mass  analysis. 

a.  August  1988  alongshore  geostrophic  velocities  (seasonal 
normal) 

For  the  two  week  period  leading  up  to  the  cruise,  the  alongshore 
component  of  surface  wind  stress  (Figure  33)  was  essentially  zero  with  a  brief 
period  of  equatorward  stress  on  July  29th.  Stations  along  the  POST  were 
occupied  before  any  significant  equatorward  component  in  the  wind  stress 
developed.  The  vertical  section  of  alongshore  geostrophic  velocity  (Figure  34) 
depicts  equator  surface  flow  between  stations  9  and  22  associated  with  the 
CC.  The  core  (v  <  -20.0  cm  s_1)  located  near  station  19  is  displaced  30  km 
inshore  of  the  salinity  minimum  (<  32.8  PSU,  Figure  14).  The  region  between 
stations  19  and  21  may  coincide  with  the  core  of  the  CC,  which  carries  low 
salinity  water  southward  and  may  occur  in  bands  defined  by  one  or  two  pairs  of 
stations  separated  by  up  to  75-150  km  (Lynn  and  Simpson  1987).  Subsurface 
poleward  flow  was  observed  predominantly  inshore  of  station  15.  The  CUC  is 
represented  by  two  cores  of  poleward  flow  located  at  stations  3  and  9,  between 
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Figure  33.  Alongshore  component  of  wind  stress  during  cruise  CUC- 

August  1988:  (a)  Wind  stress  data  from  buoy  46042,  Monterey  Bay; 
(b)  Wind  stress  data  from  buoy  46028,  Cape  San  Martin.  The  cruise 
period  is  indicated  by  the  solid  horizontal  line  (8/3  =  JD  215;  8/8  =  JD 
220). 
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Figure  34.  Vertical  section  (0-1000  dbar)  of  alongshore  geostrophic 
velocity  for  cruise  CUC-August  1988:  The  contour  interval  is 
5.0  cm  s'1.  Dashed  lines  equatorward  and  solid  lines  are  poleward. 


12  and  42  km  of  the  coast  at  a  depth  of  190  m.  Velocities  in  excess  of  20  cm  s'1 
were  observed  in  the  core  offshore,  while  inshore  they  exceeded  35  cm  s_1.  The 
position  of  these  cores  agreed  favorably  with  the  region  of  positive  spiciness 
anomaly  shown  in  Figure  16. 

Strong  poleward  flow  (v  >  60  cm  s-1)  at  the  surface  was  observed 
within  20  km  of  the  coast.  Such  inshore  countercurrents  may  be  driven  by  a 
residual  pressure  gradient,  causing  a  northw'ard  acceleration  during  wind 
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relaxations,  however  the  dynamics  involved  are  not  clearly  understood  (Huyer  et 
al.  1989).  Wickham  (1975)  observed  drifts  in  excess  of  50  cm  s-1  using 
subsurface  drouges  (at  50  m  depth)  deployed  in  August  1972  between  Point  Sur 
and  Cypress  Point,  which  he  attributed  to  non-geostrophic  components  in  the 
flow  field.  The  33.4-33.7  PSU  isohalines  surfaced  between  stations  3  and  5 
(Figure  14)  depicting  a  frontal  structure  not  present  in  the  temperature  field. 
This  developed  after  an  upwelling  event  and  surface  heating  erased  the 
temperature  signal.  It  is  not  clear  how  this  frontal  structure  of  salinity  is  related 
to  the  strong  poleward  flow  observed  between  these  stations. 

b.  May  1988  alongshore  geostrophic  velocities  (weak 

poleward  flow) 

This  cruise  was  preceded  by  a  period  of  strong  equatorward  winds, 
at  times  in  excess  of  23  knots.  The  resulting  alongshore  component  of  surface 
wind  stress  (Figure  35)  reached  or  exceeded  -2.0  dyne  cm'2  and  persisted 
through  the  first  two  days  of  the  cruise.  A  vertical  section  of  alongshore 
geostrophic  velocity  from  0-1000  dbar  (Figure  36)  depicts  equatorward  flow 
above  200  in  over  most  of  the  POST,  expected  with  strong  equatorward  winds. 
An  equatorward  surface  jet  (v  <  -25.0  cm  s'1)  is  present  between  stations  3  and  4, 
while  velocities  in  excess  of  35  cm  s’1  are  seen  inshore  at  station  1.  These 
features  are  characteristic  of  coastal  upwelling  which  produces  an  equatorward 
jet  in  response  to  strong  horizontal  density  gradients  resulting  from  the  raising 
of  isopycnals  near  the  coast  (Huyer  1983).  Upward  sloping  isopycnals  and 
horizontal  density  fronts  are  seen  in  the  vertical  section  of  density  (Figure  17) 
and  correspond  to  the  regions  of  intense  equatorward  flow  (Figure  36). 
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Figure  35.  Alongshore  component  of  wind  stress  during  the  May 
1988  Student  Cruise:  (a)  Wind  stress  data  from  buoy  46042, 
Monterey  Bay;  (b)  Wind  stress  data  from  buoy  46028,  Cape  San 
Martin.  The  cmise  period  is  indicated  by  the  solid  horizontal  line  (5/4 
=  JD  124;  5/8  =  JD  128). 
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A  secondary  subsurface  core  of  equatorward  flow  (v  <  -10.0  cm  s  ’)  was 
observed  below  station  5,  with  velocities  of  -5.0  cm  s_1  reaching  750  m.  This 
deep  equatorward  flow  is  too  deep  for  direct  wind  forcing  and  its  cause  is 
unknown. 


STM  AY  1988  (E) 
GE03  YEL  (CM/S) 


Figure  36.  Vertical  section  (0-1000  dbar)  of  alongshore  geostrophic 

velocity  for  the  May  1988  Student  Cruise:  The  contour  interval 
is  5.0  cm  s_1.  Dashed  lines  are  equatorward  and  solid  lines  are 
poleward. 

Weak  poleward  flow  was  observed  below  200  m  inshore  of  station 
4  with  a  maximum  of  4.5  cm  s*1  near  460  m.  Due  to  the  contouring  interva'  this 
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region  is  not  depicted  in  Figure  34.  The  location  of  this  poleward  flow  is 
indicative  of  the  CUC  which  historically  lies  below  the  surface  at  this  time  of 
year.  Its  anomalously  weak  magnitude  (<  5  cm  s_1)  may  result  from  the  intense 
equatorward  winds  which  preceded  this  cruise.  Poleward  flow  (>  5.0  cm  S'1) 
was  also  observed  between  stations  7  and  1 1  with  very  weak  equatorward  flow  in 
the  upper  50  m  between  station  8  and  9. 

c.  November  1988  alongshore  geostrophic  velocities  (strong 

mesoscale  feature ) 

The  water  mass  analysis  showed  that  an  anticyclonic  mesoscale  eddy 
was  the  dominant  feature  on  this  cruise,  which  is  related  to  the  geostrophU 
velocities  (Figure  37).  Poleward  velocities  in  excess  of  55  cm  s_1  occur  at 
station  1 1  near  1 10  m  depth.  A  second  subsurface  maximum  occurs  at  station  13 
w'ith  velocities  in  excess  of  40  cm  s*1. 

Inshore  of  station  5,  the  CUC  was  observed  at  the  surface  and  along 
the  continental  slope,  with  a  core  depth  of  460  m.  It  appears  the  mesoscale  eddy 
has  limited  the  horizontal  extent  of  the  undercurrent  to  within  10  km  of  the 
slope.  The  loct  n  of  the  CUC  in  Figure  37  coincides  with  the  downward 
sloping  isopycnals  shown  in  Figure  19  which  result  from  the  warmer,  more 
saline  waters  of  the  CUC.  The  positive  values  of  spiciness  anomaly  which  occur 
within  the  mesoscale  feature  are  likely  the  result  of  lateral  entrainment  of  CUC 
water  as  discussed  earlier. 

Poleward  flow'  was  observed  from  station  16  to  the  midpoint 
between  stations  18  and  19.  From  earlier  analysis  of  the  water  mass  structure,  it 
was  found  that  throughout  this  region  the  water  was  of  subarctic  origin.  It  is 
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Figure  37.  Vertical  section  (0-1000  dbar)  of  alongshore  geostrophic 

velocity  for  cruise  CUC-November  1988:  The  contour  Interval 
is  5.0  cm  S'1.  Dashed  lines  are  equatorward  and  solid  lines  are 
poleward. 

possible  that  wind  relaxations  prior  to  and  during  the  cruise  (Figure  38)  was  the 
mechanism  b  .hind  this  poleward  flow,  and  poleward  advection  of  subarctic 
water,  normally  associated  with  equatorward  flow.  Offshore  of  this  poleward 
flow  was  the  equatorward  flow-ing  CC,  with  velocities  less  then  -20.0  cm  s1. 
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d.  February  1989  alongshore  geostrophic  velocities  ( strong 
equatorward  flow  near  coast) 

Wind  data  for  this  cruise  could  only  be  obtained  from  the  Cape  San 

* 

Martin  buoy,  since  the  Monterey  Bay  buoy  was  inoperative  during  the  last  half 
of  January  and  all  of  February.  The  record  from  buoy  46028  (Figure  39)  shows 
two  strong  equatorward  wind  events  preceding  the  cruise  by  several  days,  and  a 
weaker  one  just  prior  to  the  start  of  the  cruise.  A  relaxation  period  of  little  cr 
no  wind  occurred  four  days  before  the  start  of  the  cruise  and  another  at  the  start 
of  the  cruise.  By  February  4th  the  winds  were  once  again  equatorward  for  two 
days,  switching  to  poleward  as  time  progressed. 

The  vertical  section  of  alongshore  geostrophic  velocity  (Figure  40) 
shows  the  CUC  at  the  surface  between  station  5  and  1 1  extending  to  a  depth  of 
600  m.  It  has  a  subsurface  core  (v  >  35.0  cm  s'1)  lying  near  a  depth  of  100  m, 
15  km  from  the  coast.  Offshore  of  the  CUC,  the  flow  is  generally  equatorward 
with  large  areas  of  weak  flow  between  0.0  and  -5.0  cm  s'1.  Correlation  of  wind 
records  between  buoys  46028  and  46042  for  each  cruise  show  stronger 
equatorward  wind  events  at  buoy  46028  to  the  south  and  stronger  poleward  wind 
events  at  buoy  46042  to  the  north.  If  it  is  assumed  a  similar  pattern  existed  on 
this  cruise,  the  equatorward  wind  event  just  prior  to  the  cruise  may  have  been 
stronger  than  represented  in  Figure  39.  The  poleward  flow  observed  between 
stations  5  and  1 1  would  then  be  the  result  of  the  relaxation  of  the  equatorward 
winds.  This  may  also  account  for  the  stronger  flow  of  the  CUC. 

Inshore  of  station  5  is  a  region  of  strong  equatorward  flow  (v  < 
-40.0  cm  s_1J.  A  salinity  minimum  (S  <  33.3  PSU)  existed  between  stations  5  and 
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Figure  38.  Alongshore  component  of  wind  stress  during  cruise  CUC- 
November  1988:  (a)  Wind  stress  data  from  buoy  46042,  Monterey 
Bay;  (b)  Wind  stress  data  from  buoy  46028,  Cape  San  Martin.  'Die 
cruise  period  is  indicated  by  the  solid  horizontal  line  (1 1/14  =  JD  318; 
11/19  =  JD  323). 
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Figure  39.  Alongshore  component  of  wind  stress  during  cruise  CUC- 
February  1989:  Wind  stress  data  from  buoy  46028,  Cape  San 
Martin.  The  cruise  period  is  indicated  by  the  solid  horizontal  line  (2/3  = 
JD  34;  in  =  JD  38). 

6  corresponding  to  the  density  minimum  shown  in  Figure  41.  The  upward 
sloping  isopycnals  near  the  coast  suggest  that  the  salinity  minimum  may  be  the 
mechanism  behind  the  strong  equatorward  flow  inshore  of  the  CUC.  The 
salinity  minimum  may  be  the  result  of  ambient  surface  water  trapped  between 
previously  upwelled  water  and  higher  salinity  water  near  the  coast. 

e.  May  1989  alongshore  geostrophic  velocities  (poleward 
flow  near  coast) 

The  alongshore  component  of  wind  stress  was  between  -1.0  to  -2.0 
dyne  enr2  three  days  prior  to  this  cruise  and  continued  through  the  first  two  days 
of  the  cruise  (Figure  42).  A  pattern  of  poleward  and  equatorward  flows  is 
depicted  in  the  alongshore  geostrophic  velocities  (Figure  43)  .  The  most  intense 
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Figure  40.  Vertical  section  (0-1000  dbar)  of  alongshore  geostrophic 

velocity  for  cruise  CUC-February  1989:  The  contour  interval  is 
5.0  cm  s-'.  Dashed  lines  are  equatorward  and  solid  lines  are  poleward. 
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Figure  42.  Alongshore  component  of  wind  stress  during  the  May 
1989  Student  Cruise:  (a)  Wind  stress  data  from  buoy  46042, 
Monterey  Bay:  (b)  Wind  stress  data  from  buoy  46028,  Cape  San 
Martin.  The  cruise  period  is  indicated  by  the  solid  horizontal  line  (5/4 
=  JD  124;  5/8  =  JD  128). 
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Figure  43.  Vertical  section  (0-1000  dbar)  of  alongshore  geostrophic 
. elocity  for  the  May  1989  student  cruise:  The  contour  interval 
is  5.0  cm  s‘‘.  Dashed  lines  are  equatorward  and  solid  lines  are 
poleward. 

flows  occur  where  surface  density  fronts  are  present,  stations  7  through  1 1  and 
14  through  16.  A  surface  equatorward  jet  (core  v  <  -40.0  cm  s_1,  v  <  -20.0  cm 
S'1  at  150  m)  is  coincident  with  the  density  front  between  stations  7  through  1 1 
(Figure  26).  The  deeper  sliucturc  penetrating  to  a  depth  of  approximately  450 
m  coincides  with  the  depression  of  the  isopycnals  below  200  m  (Figure  26).  The 
surface  jet  is  a  characteristic  feature  of  an  upwclling  regime  (Huycr  1983).  An 
equatorward  (v  <  -15.0  cm  s1)  jet  and  weaker  poleward  jet  (v  >  5.0  cm  s_1) 


coincide  with  the  frontal  zone  between  station  14  and  16.  This  entire  region  is 
believed  to  be  associated  with  the  CC,  with  poleward  flow  resulting  from  the 
isopycnal  slopes  (downward  toward  the  coast). 

Inshore  of  station  5  is  a  region  of  poleward  flow  which  extends 
from  the  surface  to  a  depth  of  750  m  along  the  continental  slope.  In  the  presence 
of  upwelling  favorable  winds  the  CUC  is  expected  to  lie  at  depth  along  the 
continental  slope  with  no  surface  signature  (Chelton  1984).  However,  due  to  the 
presence  of  a  salinity  minimum  near  the  surface  (Figure  24)  between  stations  3 
and  5,  poleward  flow  extends  to  the  surface.  Based  upon  the  distribution  of 
salinity,  the  CUC  is  believed  to  lie  beneath  and  is  dynamically  distinct  from  the 
poleward  flow  at  the  surface.  Isohalines  at  depth  slope  upward  toward  the 
continental  slope  characteristic  of  the  more  saline  waters  of  the  undercurrent, 
while  those  at  the  surface  illustrate  the  presence  of  a  salinity  minimum.  Surface 
temperatures  between  station  4  and  6  were  greater  than  13°C,  making  this 
intrusion  warmer  and  less  saline  than  surrounding  waters.  This  warm  region  is 
also  reflected  in  SST  from  NOAA  1 1  AVHRR  satellite  imagery  (Figure  44). 
The  uppe"  -  ater  column  undergo00  e°ostrophic  adjustment  resulting  in  near 
surface  poleward  flow.  The  low  salinity  water  may  well  be  ambient  surface 
water  trapped  between  freshly  upw'elled  water  inshore  and  a  previously  existing 
front  offshore  from  the  equatorward  ''dnd  event  and  relaxation  six  days  prior  to 
the  cruise. 

/.  July  1989  alongshore  geostrophic  velocities  (anomalous 
poleward  flow) 

Winds  for  the  14  day  period  leading  up  to  this  cruise  were  predominantly 
equatorward  with  a  substantial  wind  relaxation  lasting  approximately  three  days 
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just  two  days  before  the  cruise.  Winds  during  the  cruise  were  equatorw'ard  at 
both  buoy  locations  with  alongshore  wind  stress  between  -0.75  dyne  cm  2  at  buoy 
46042  and  -2.5  dyne  cnr2  at  buoy  46028  (Figure  45).  The  vertical  section  of 
alongshore  geostrophic  velocity  (Figure  46)  is  dominated  by  poleward  flow 
along  most  of  the  transect.  Weak  equatorward  flow  (v  >  -5.0  cm  s_1)  is  found  at 
the  surface  offshore  of  station  17  and  the  -5.0  cm  s'1  contour  becomes  visible 
near  the  surface  at  station  21.  There  is  weak  equatorward  flow  between  stations 
1  and  2  and  8  and  9  in  response  to  equatorward  winds  however  it  is  minor 
compared  to  the  poleward  flow. 

The  strong  poleward  flow  throughout  the  section  is  believed  to 
result  from  the  substantial  wind  relaxation  which  preceded  this  cruise.  When 
equatorward  winds  relax,  a  residual  alongshoie  pressure  gradient  can  produce  a 
northward  acceleration.  The  wind  records  indicate  that  the  relaxation  at  the 
southern  buoy  was  of  longer  duration  than  to  the  north.  In  their  study  of 
variability  along  the  California  shelf  induced  by  local  and  remote  winds  during 
the  Coastal  Ocean  Dynamics  Experiment  (CODE),  Davis  and  Bogden  (1989) 
discuss  a  length  scale  500  km  associated  with  remote  forcing.  Without  wind 
records  from  bu'»v  i, other  south  it  can  only  be  hypothesized  that  this  poleward 
flow'  was  in  response  to  a  large  scale  wind  relaxation  which  might  have  occurred 
previously  over  the  entire  central  and  southern  California  region.  The  CUC  was 
observed  to  occupy  the  region  of  the  continental  slope  between  200  m  600 
m,  with  a  core  depth  of  400  m.  Above  the  CUC.  poleward  flow'  is  believed  to  be 
associated  with  the  wind  relaxation  discussed  earlier.  The  weaker  equatorward 
flow  offshore  may  be  the  oceanic  response  to  the  dominant  equatorward  winds. 


on 


Figure  44.  NOAA  AVIIKR  satellite  image  from  2044  (IMF, 

3  May  1989:  Darker  shades  near  the  coast  depict  colder,  more 
recently  upwellcd  waters.  Lighter  shades  found  offshore  depict 
warmer  water.  SS  I  increases  toward  the  offshore  corner  of  the 
imace.  Grid  spacine  equals  30'.  Note  the  warm  region  between 
122225°- 123.00°  \\\  and  35.75°-36.50  ’  N. 
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Figure  45.  Alongshore  component  of  wind  stress  during  cruise  CUC-July 
19X9:  (a)  Wind  stress  data  from  buoy  46042,  Monterey  Bay; 

(h)  W  ind  stress  data  from  buoy  46028,  Cape  San  Marin.  The  cruise 
period  is  indicated  bv  the  solid  horizontal  line  (7/28  =  JD  209;  8/3  = 

JI)  215). 
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Figure  46.  Vertical  section  (0-1000  dbar)  of  alongshore  gcostrophic 
velocity  for  cruise  CUC-JuIy  1989:  The  contour  interval  is 
5.0  cm  s'1.  Dashed  lines  are  equatorward  and  solid  lines  are  polew  ard. 

This  would  explain  why  the  poleward  flow  becomes  both  deeper  and  weaker 
offshore  as  flow  conditions  begin  to  return  to  normal.  Comparisons  with  similar 
observations  from  July  1981  and  1984  (Chelton  et  al.  1988)  follow  in  Chapter  V. 
g.  November  1989  alongshore  gcostrophic  velocities  (deep 
equatorward  flow ) 

The  alongshore  component  of  surface  stress  (Figure  47)  shows  a 
series  of  equatorward  wind  events  and  relaxations  prior  to  and  during  this 
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cruise.  In  opposition  to  the  earlier  cruises  where  the  equatorward  winds  were 
stronger  at  the  southern  buoy,  the  Monterey  Bay  buoy  reflects  greater  values  of 
equatorward  stress  than  the  Cape  San  Martin  buoy.  Unlike  the  previous  six 
cruises,  the  chc;ce  of  1000  dbar  as  a  level  of  no  motion  appears  to  be  in  error  on 
this  cruise.  Focussing  on  the  upper  1000  m  (Figure  48),  the  CC,  normally  found 
within  the  upper  300  m  (Lynn  and  Simpson  1987),  has  considerable  vertical 
extent.  Inshore  of  station  19  the  flow  pattern  becomes  more  complicated. 
The  complex  pattern  shown  in  Figure  48  may  be  in  response  to  the  rapidly 
changing  wind  field  during  the  cruise.  From  water  mass  analysis,  there  was 
significant  intrusion  of  either  North  Pacific  Central  water  at  depth  or  an 
unusually  strong  recirculation  of  CUC  water  between  stations  20  and  22,  which 
forced  a  large  area  of  equatorward  flow-  (Figure  48)  extending  to  depths  of  more 
than  700  m  near  stations  20-22,  and  poleward  flow  greater  than  5.0  cm  s*1 
extending  from  1800  m  to  the  bottom  (Figure  49).  Deep  poleward  flow  of  this 
magnitude  is  unlikely  and  suggests  that  on  this  cruise  2000  dbar  may  have  been  a 
better  choice  for  the  LNM. 

The  CUC  is  observed  at  the  surface  between  stations  4  and  10,  w  ith 
a  subsurface  maximum  (v  >35  cm  s_1)  at  a  depth  of  70  m.  Inshore  of  this  is  a 
region  of  weaker  equatorward  flow  within  15  km  of  the  coast  believed  to  result 
from  the  equatorward  w-inds  present  at  the  start  of  the  cruise.  The  poleward 
flow  and  strong  horizontal  velocity  gradient  seen  near  station  18  and  19 
correspond  to  the  isopycnal  slopes  (downward  toward  the  coast)  and  surface 
front  seen  in  Figure  29. 
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Figure  47.  Alongshore  component  of  wind  stress  during  cruise  CUC- 
November  1989:  (a)  Wind  stress  data  from  buoy  46042,  Monterey 
Bay;  (b)  Wind  stress  data  from  buoy  46028,  Cape  San  Martin.  The 
cruise  period  is  indicated  by  the  solid  horizontal  line  (1 1/15  =  JD  320; 

1 1/22  =  JD  327). 
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CUC- NOVEMBER  1939  (E)  1 — 

GEOS.  VEL.  (CM/S)  10  KM 


Figure  48.  Vertical  section  (0-1000  dbar)  of  alongshore  geostrophic 

velocity  for  cruise  CUC-November  1989:  The  contour  interval 
is  5.0  cm  sb  Dashed  lines  are  equatorward  and  solid  lines  are 
poleward. 

Another  interesting  feature  found  in  all  the  sections  (Figure  49,  for 
example)  is  the  deep  banding  of  the  velocity  structure  which  occurs  below  the 
reference  level  between  stations  9  and  18.  This  banding  was  first  noticed  by 
Berryman  (1989)  in  his  study  of  currents  off  Point  Sur  during  February  1989. 
The  topography  in  this  area  in  complicated  by  the  presence  of  the  southern 
extension  of  the  Monterey  canyon  which  is  responsible  for  the  steepness  of  the 
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continental  slope  in  this  area.  This  banded  structure  is  believed  to  by  the  result 
of  high  frequency  motions  interacting  with  the  steep  canyon  topography.  The 
flow  is  not  in  geostrophic  balance  as  the  magnitude  of  flow  sometimes  reached  or 
exceeded  40  cm  s'1  (November  1988,  not  shown).  The  magnitudes  of  these  flows 
are  far  greater  than  those  resulting  from  error  in  the  calculations  supporting  the 
idea  that  a  high  frequency  motion  ex'  ,s,  disturbing  the  mass  field  to  produce  this 
banded  structure,  which  requires  further  investigation. 

For  the  sake  of  uniformity,  the  reference  level  was  lett  ai  1000 
dbar  for  comparison  with  earlier  cruises.  By  using  a  deeper  level  one  would 
expect  to  see  an  increase  in  the  magnitude  of  the  values  shown  in  Figure  48  and 
the  elimination  of  the  deep  poleward  flow  seen  in  Figure  49. 
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Figure  49.  Vertical  section  (0-4500  dbar)  of  alongshore  geostrophic 

velocity  for  cruise  CUC-Novembcr  1989:  The  contour  interval 
is  5.0  cm  s*1.  Dashed  lines  are  equatorward  and  solid  lines  are 
poleward. 
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V.  DISCUSSION 


This  Chapter  focuses  on  similarities  and  differences  between  the  observations 
of  this  study  and  those  of  earlier  studies  conducted  off  Point  Sur,  California, 
with  emphasis  placed  on  the  possible  dynamics  of  the  flow.  As  pieviously 
mentioned,  most  of  the  earlier  studies  off  Point  Sur  utilized  many  years  of 
hydrographic  data  collected  in  connection  with  the  CalCOFI  program,  and  our 
observed  differences  may  in  most  cases  lie  within  the  standard  deviations  from 
their  seasonal  averages. 

A.  WATER  MASS  CHARACTERISTICS 

The  water  masses  present  off  Point  Sur  as  indicated  by  the  spiciness  anomaly 
appear  to  be  consistent  with  those  discussed  by  Lynn  and  Simpson  (1987)  based 
on  the  CalCOFI  data  set.  The  core  of  the  CC  is  characterized  by  relatively  low 
temperature  and  low  salinity.  The  western  boundary  of  the  CC  as  defined  by 
Lynn  and  Simpson  (1987)  is  a  salinity  frontal  structure  which  lies  between  the 
southern  extension  of  the  Subarctic  Frontal  Zone  and  waters  of  the  eastern  North 
Pacific.  Along  CalCOFI  line  70  (Figure  50),  the  mean  salinity  for  July  along  the 
surface  where  ot  =  25.0  shows  a  minimum  of  less  than  33.0  PSU  between  124° 
and  125°  W.  Observations  from  May  1988  and  November  1988  (not  shown)  also 
reveal  salinities  less  33.0  PSU  near  the  surface  in  the  CC,  while  during  May  1989 
(Figure  24)  and  November  1989  (not  shown),  values  less  than  32.9  PSU  were 
observed.  Salinity  values  less  than  32.8  PSU  were  observed  in  the  CC  during 
August  1988  (Figure  14).  The  Coastal  Transition  Zone  (CTZ)  sampling  grid 
covers  a  section  of  the  northern  California  coast  approximately  200  km  in 
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length,  including  Point  Arena  and  Point  Reyes,  and  extends  to  about  200  km 
offshore  (Huyer  et  al.  1990).  Surface  salinities  found  in  the  offshore  region  of 
the  CTZ  sampling  grid  in  the  summer  of  1988  were  commonly  less  than  33.0 

t 

PSU  with  a  low  salinity  core  (S  <  32.6  PSU)  extending  offshore  from  Point 
Arena  coincident  with  a  strong  baroclinic  jet  (Huyer  et  al.  1990).  During  the 
survey  period  a  sharp  meander  developed  in  this  jet.  The  lower  salinity  values 
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Figure  50.  Mean  salinity  on  the  surface  where  Ct  =  25.0  for  July: 

The  data  presented  here  comes  from  the  CalCOFl  data  set. 

The  bold  line  represents  the  intersection  of  this  density  surface 
and  the  sea  surface.  The  density  is  less  inshore  of  the  bold  line. 
Contour  interval  is  0. 1  PSU.  (Source:  Lynn  and  Simpson  198  7) 
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observed  along  the  POST  are  characteristic  of  waters  farther  to  the  north  and 

may  result  from  the  southern  transport  of  fresher  water  through  energetic  jets 

and  meanders  present  within  the  transition  zone  as  described  by  Lynn  and 

0 

Simpson  (1987)  and  the  Coastal  Transition  Zone  Experiment. 

Unusually  high  surface  salinities  (>33.5  PSU)  were  observed  during  the  July 
1989  cruise,  nearly  0.6  PSU  higher  than  the  August  1988  cruise.  Chelton  et  al. 
(1988)  also  found  anomalously  high  salinity  values  off  Point  Sur  in  July  1981 
and  July  1984.  Salinities  at  the  10  m  level  for  both  these  periods  (Figure  51) 
show  intrusions  of  higher  salinity  water  from  the  south  in  opposition  to  the 
surface  equatorward  winds,  suggestive  of  remote  forcing  from  the  south.  The 
vertical  sections  of  salinity  from  August  1988  (Figure  14)  and  July  1989  (Figure 
Z6),  illustrate  that  the  surface  values  were  much  higher  in  1989  and  appear  to 
mimic  the  conditions  observed  by  Chelton  et  al.  (1988). 

The  deep  intrusion  of  warm,  salty  water  and  accompanying  deep  equatorward 
flow  in  November  1989  is  a  significant  deviation  from  the  normal  seasonal 
picture  of  this  area.  Below  200  m  over  the  entire  POST  the  isopycnals  (Figure 
29)  have  a  domed  appearance  centered  at  station  19  and  slope  downward  both 
offshore  and  inshore.  This  doming  of  the  isopycnal  surfaces  is  common  during 
winter  off  the  California  coast  (Reid  1973).  The  intrusion  of  the  warmer,  more 
saline  water  may  represent  either  North  Pacific  Central  water  which  usually  lies 
farther  offshore  (128°  W  at  this  latitude,  Lynn  and  Simpson  1987)  or  an 
unusually  strong  recirculation  of  CUC  water.  Similar  but  weaker  intrusions 
were  also  observed  in  August  1988,  May  1989,  and  July  1989  and  may  again 
result  from  either  of  these  sources.  To  determine  the  exact  cause  a  more 
intensive  T-S  analysis  is  required. 
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Figure  51.  Horizontal  sections  of  salinity  at  10  m  depth  from  CalCOFI 
data:  The  top  section  illustrates  the  July  seasonal  average;  the  center 
section  illustrates  the  conditions  for  July  1981;  the  bottom  section 
illustrates  the  conditions  for  July  1984.  (Source:  Chelton  et  al.  1988) 

B.  ALONGSHORE  GEOSTROPHIC  VELOCITIES 

In  his  study  of  alongshore  geostrophic  velocities  along  the  central  California 
coast,  Chelton  (1984)  used  23  years  of  hydrographic  data  along  CalCOFI  lines  70 


and  80,  off  Point  Sur  and  Point  Conception,  respectively,  to  establish  the  seasonal 
variability  of  alongshore  geostrophic  currents  across  these  sections.  His  analysis 
was  restricted  to  stations  where  the  bottom  depth  exceeded  500  m,  which  was 
roughly  within  15  km  of  the  coast  off  Point  Sur.  These  historical  data  (Figure 
52)  reveal  the  CC  within  the  upper  200  m  flowing  equatorward  with  its  core 
located  between  100  and  200  km  off  the  coast.  There  are  two  maxima  during  the 
year,  between  February  and  March  and  July  through  September  where  the 
velocities  exceed  -9.0  cm  s'1.  A  narrow  second  maximum  occurs  near  the  coast 
during  March-April  and  July-September  with  velocities  exceeding  -7.5  cm  s'1 
and  -5.0  cm  s1,  respectively.  Chelton  (1984)  indicates  this  narrow  jet  may  be 
even  more  intense  over  the  continental  shelf.  This  banded  flow  structure  was 
also  observed  by  Lynn  and  Simpson  (1987),  who  found  that  equatorward  flow  in 
excess  of  -4  cm  s_1  rarely  occurred.  The  CUC  off  Point  Sur  (Figure  52)  appears 
to  be  confined  to  the  continental  slope  within  75-100  km  of  the  coast,  with  no 
nearshore  poleward  flow  at  depth  in  March  through  May.  Poleward  flow 
extends  to  the  surface  from  October  through  February  (the  Davidson  Current?) 
with  a  maximum  surface  velocity  of  14.0  cm  s1  in  December,  and  is  subsurface 
throughout  the  remainder  of  the  year.  Chelton  points  out  that  while  seasonal 
poleward  flow  is  small  near  500  dbar,  it  may  not  necessarily  be  zero  suggesting 
the  use  of  500  dbar  as  the  LNM  may  be  in  error.  Based  upon  the  observations  of 
our  study  the  use  of  500  dbar  as  the  LNM  is  clearly  in  error.  Equatorward  flow 
associated  with  the  CC  was  observed  below  500  dbar  on  several  occasions  and  the 
nearshore  poleward  flow  of  CUC  was  found  as  deep  as  600  m  in  some  cases. 

Peak  velocities  of  the  CC  and  CUC  along  the  POST  typically  exceeded  those 
observed  during  CalCOFI  (Figure  52).  The  CC  had  velocities  in  excess  of  -20.0 
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cm  s'!  in  May  1988,  November  1988  and  May  1989.  The  CUC  had  velocities  in 
excess  of  35  cm  s_1  in  August  1988,  February  1989  and  November  1989.  These 


Figure  52.  Vertical  section  (0-400  m)  of  seasonal  alongshore  geostrophic 
velocity  relative  to  500  dbar  along  CalCOFI  line  70  off  Point 
Sur:  Shaded  regions  depict  poleward  flow.  The  contour  interval  is 
2. 5  cm  s' 1.  (Source:  Chelton  1984) 

higher  geostrophic  velocities  were  likely  the  result  of:  1)  the  decreased  station 
spacing  of  the  POST;  and  2)  the  deeper  LNM  used  in  this  study.  Chelton  (1984) 
points  out  that  the  coarse  station  spacing  associated  with  the  CalCOFI  sampling 
grid  tends  to  provide  a  smoothed  version  of  the  flow  field  and  potentially  misses 
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narrow,  higher  velocity  jets  near  the  coast.  The  decreased  station  spacing 
inshore  provides  better  definition  of  these  narrow  jets  and  increases  the 
magnitude  of  the  mean  flow.  Significant  vertical  shear  stilj  existed  at  and  below 
500  dbar,  which  suggests  that  the  deeper  LNM  contributed  to  the  larger  surface 
velocities  observed  in  this  study. 

Based  upon  the  analysis  of  water  mass  characteristics  and  geostrophic 
velocities,  the  spatial  extent  and  peak  core  velocities  of  the  CUC  for  each  of  the 
cruises  have  been  defined  (Table  8).  The  core  of  the  undercurrent  was  observed 
within  42  km  of  the  coast  and  between  70  m  to  460  m  depths  on  all  cruises, 
which  generally  agrees  with  the  observations  of  Wickham  et  al.  (1987).  These 
values  are  much  greater  than  those  shown  in  Figure  50  and  may  be  the  result  of 
decreased  station  spacing  near  the  coast.  The  CUC  was  present  at  the  surface  in 
February  1989  and  November  1989  and  subsurface  throughout  the  remainder  of 
the  year,  similar  to  the  results  of  Chelton  (1984).  Poleward  flow  was  only 
present  at  the  surface  inshore  of  station  3  (Figure  37)  due  to  the  presence  of  a 
strong  mesoscale  feature  in  November  1988,  and  was  below  this  feature 
offshore.  Wickham  et  al.  (1987)  observed  increased  temperatures  and  salinities 
nearshore  in  winter  within  the  upper  layers.  The  warmest  and  saltiest  waters 
were  found  at  the  depth  of  the  CUC  in  agreement  with  the  observations  of 
W'ooster  and  Jones  (1970)  off  northern  Baja  California,  who  also  observed  the 
CUC  at  depth  below  shallow  weak  equatorward  and  sometimes  poleward  flow  in 
summer.  From  their  observations  off  northern  Baja  near  3l°N,  Wooster  and 
Jones  (1970)  found  the  CUC  to  lie  close  to  continental  slope.  It  had  a  width  of 
approximately  20  km,  a  thickness  of  300  m,  and  velocities  approaching 
30  cm  S'1.  These  earlier  studies  agree  tavoraWy  with  the  observations  of  this 
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study  although  their  measurements  were  taken  off  Cape  San  Martin,  roughly  70 
km  south  of  Point  Sur,  and  northern  Baja,  roughly  800  km  to  the  south. 

Table  8.  SPATIAL  EXTENT  AND  CORE  VELOCITIES  OF  THE 
CALIFORNIA  UNDERCURRENT:  Based  upon  the 
observations  from  seven  cruises  along  the  Point  Sur  Transect 
between  May  1988  and  November  1989. 


Cruise 

Distance  from 
shore  (km) 

Core  depth  (m) 

Maximum 
Velocity  (cm  s'1) 

STMAY  1988 

22 

460 

<  5 

CUC-Aucust  1988 

12/42 

190 

>  35  / >  20 

CUC-November  1988 

20 

460 

>  25 

CUC- February  1989 

15 

100 

>  35 

STMAY  1989 

17 

160 

>20 

CUC-July  1989 

23 

400 

>  10 

CUC-November  1989 

28 

70 

>35 

indicating  properties  of  the  CUC  are  not  observed  to  vary  widely  along  the 
California  coast. 

The  strong  \  <eward  surface  flow  observed  in  May  1989  inshore  of  station  5 
(Figure  43)  appeared  to  be  the  result  of  warmer,  fresher  water  at  the  surface, 
not  associated  with  the  CUC,  which  was  observed  at  depth.  This  poleward  flow' 
resulted  from  the  geostrophic  adjustment  of  the  water  column  in  response  to  the 
warm,  fresh  water  found  near  the  surface  farther  offshore.  Strong  poleward 
surface  flow  was  also  observed  in  August  1988  and  July  1989  although  the 
dynamics  involved  differ  from  the  flow  observed  in  May  1989.  In  July  1989.  a 
significant  relaxation  of  equatorward  winds  (Figure  45)  occurred  tw'o  days  prior 
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to  this  cruise  and  lasted  three  days.  The  magnitude  of  this  relaxation  was  greater 
at  the  southern  buoy,  46028. 

The  reversal  of  near  surface  currents  in  response  to  relaxation  of 
equatorward  winds  has  been  observea  using  moored  current  meters  during 
CODE  (Davis  1985).  In  May  1981,  equatorward  winds  of  10  m  s1  at  central 
CODE  stations  relaxed  to  5  m  s*1.  Within  one  day  after  this  relaxation  current 
meters  located  125  km  to  the  south  (in  the  southern  section  of  the  CODE  region) 
exhibited  a  reversal  in  flow  from  equatorward  to  poleward.  Three  days  later 
this  flow  reversal  was  present  in  current  meters  70  km  to  the  north  of  the 
southern  current  meters,  indicating  poleward  propagation  of  this  event.  Davis 
and  Bogden  (1989)  discuss  the  effect  of  local  and  remote  wind  forcing  on  shelf 
circulation  in  the  CODE  region.  In  their  study,  coastally  trapped  waves  depicted 
in  sea  level  records  were  found  to  be  coherent  with  coastal  winds  500  km  to  the 
south.  It  is  hypothesized  that  the  strong  poleward  flow  exhibited  in  the  July  1989 
alongshore  geostrophic  velocities  (Figure  46)  was  in  response  to  remote  forcing 
from  the  south,  and  the  weak  equatorward  flow  and  deeper  poleward  flow 
offshore  represent  the  oceanic  response  to  the  onset  of  equatorward  winds. 
Further  evidence  to  support  this  hypothesis  comes  from  Chelton  et  al.  (1988)  as 
discussed  earlier.  Salinities  at  the  10  m  level  for  July  1981  and  July  1984  were 
higher  than  those  from  the  seasonal  average  for  July  from  the  CalCOFI  data  set 
and  are  characteristic  of  those  found  fanher  south,  indicative  of  poleward  flow. 

In  August  1988  a  region  of  intense  poleward  flow  was  observed  at  the 
surface  within  20  km  of  the  coast  (Figure  34).  A  brief  period  of  intense 
equatorward  winds  preceded  the  cruise  by  four  days  (Figure  33),  followed  by  a 
calm  period  lasting  two  to  three  days  into  the  cruise.  Present  understanding 
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suggests  that  inshore  countercurrents  may  in  part  be  driven  by:  1)  the  wind  stress 
curl  and  2)  alongshore  pressure  gradients  which  are  set  up  to  balance  strong 
equatorward  winds.  When  the  winds  relax,  a  residual  pressure  gradient 
produces  a  northward  acceleration  and  hence  poleward  flow  (Huyer  et  al.  1989). 
Drogue  trajectories  from  August  1972  depict  strong  poleward  (v  >  50  cm  s4) 
flow  inshore  between  Point  Sur  and  Cypress  Point  which  was  much  stronger  than 
the  corresponding  geostrophic  currents.  Due  to  the  design  of  these  early  drogues, 
they  were  not  truly  Lagrangian  (move  along  with  currents),  but  were  subject  to 
wind  forcing  as  well,  which  would  produce  the  higher  speeds  observed. 
Additionally,  deformations  of  the  reference  level  due  to  internal  waves  may  also 
contribute  to  the  observed  difference  between  the  flow  observed  by  drouges  and 
geostrophy  (Wickham  1975).  A  possible  explanation  for  the  unusually  intense 
poleward  flow  observed  in  our  study  is  that  it  may  be  a  response  to  the 
relaxation  of  equatorward  winds,  however  the  relaxation  occurred  four  days 
prior  to  the  cruise,  a  time  scale  longer  than  expected  from  earlier  studies.  The 
exact  cause  is  not  known. 

The  intense  equatorward  flow  found  inshore  of  station  5  during  the  February 
1989  cruise  (Figure  40)  requires  further  discussion.  Velocities  in  excess  of  -40 
cm  s4  were  observed  in  a  time  of  weak  or  no  winds  (Figure  39).  The  vertical 
section  of  salinity  depicts  a  salinity  minimum  (S  <  33.3  PSU)  which  corresponds 
to  the  minimum  observed  in  the  density  Field  (Figure  41)  between  stations  4  and 
7.  Salinity  minima  in  coastal  regions  during  winter  usually  result  from  excess 
river  run-off,  however  drought  conditions  have  existed  throughout  central  and 
southern  California  since  1988.  Therefore,  this  lower  salinity  water  is  not 
expected  to  be  the  result  of  river  run-off.  An  equatorward  wind  event  which 


occurred  just  prior  to  the  cruise  and  subsequent  offshore  transport  may  be 
responsible  for  the  weak  frontal  structure  seen  in  the  density  field.  It  is 
hypothesized  that  the  lower  salinity  water  was  just  ambient  surface  water  trapped 
between  previously  upwelled  water  and  higher  salinity  water  found  near  the 
coast.  The  observed  equatorward  flow  would  then  have  resulted  from  the 
geostrophic  adjustment  of  the  water  column  inshore.  However,  the  lack  of  wind 
data  from  the  Monterey  Bay  buoy  and  the  exact  cause  of  the  lower  salinity  water 
leave  this  hypothesis  unsubstantiated. 

The  deep  equatorward  flow  found  offshore  during  November  1989  suggests 
an  error  in  the  use  of  1000  dbar  as  the  LNM.  Due  to  the  intrusion  of  North 
Pacific  Central  water,  the  equatorward  flow  of  the  CC  penetrated  deep  within  the 
water  column.  The  deep  (ageostrophic)  flow  exhibited  near  the  steep  topography 
of  the  Monterey  Canyon  was  a  feature  present  in  all  the  cruises.  Shepard  (1975) 
indicates  that  alternating  up-  and  downcanyon  currents  with  velocities  reaching 
50  cm  s*1  can  be  found  in  submarine  canyons.  Internal  waves,  mostly  near  tidal 
frequencies,  advance  up  submarine  canyons  and  rarely  downward  (Shepard 
1975).  Strong  downcanyon  flows  may  result  from  unusual  canyon  bathymetry. 
It  is  believed  that  these  deep  ageostrophic  flows  are  the  result  of  high  frequency 
motions  interacting  with  the  steep  bathymetry  of  the  Monterey  canyon  where  it 
crosses  the  POST,  and  deserves  closer  study  in  the  future. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  purpose  of  this  study  was  to  determine  the  seasonal  variability  of 
geostrophic  velocities  and  water  mass  characteristics  along  the  Point  Sur 
Transect  (POST).  Observations  presented  throughout  this  study,  however, 
appear  to  have  revealed  interannual  rather  than  seasonal  variability.  This  is  not 
surprising  since  only  seven  cruises  were  selected  for  study  whereas  earlier 
studies  utilized  many  years  of  data  collected  along  the  CalCOFI  sampling  grid  to 
determine  the  seasonal  means.  This  study  excelled  over  the  continental  shelf  and 
slope  where  the  station  spacing  of  the  POST  is  considerably  closer  than  the 
CalCOFI  scheme.  This  decreased  spacing  enabled  the  study  of  narrow  coastal 
jets  which  could  not  be  studied  effectively  using  the  CalCOFI  scheme  (Chelton 
1984). 

On  all  cruises  the  flow  pattern  inshore  depicted  narrow  bands  of  both 
poleward  and  equatorward  flow  demonstrating  the  complexity  of  the  shelf/slope 
region  off  Point  Sur.  In  general,  observations  of  large  scale  features  such  as  the 
CC  and  in  some  cases  the  CUC  were  in  agreement  with  the  studies  of  Chelton 
(1984)  and  Lynn  and  Simpson  (1987).  The  nature  of  the  alongshore  geostrophic 
velocities  and  the  location  and  spatial  extent  of  the  undercurrent  appear  to  be 
strongly  related  to  specific  wind  events,  both  local  and  remote.  Observed 
differences  represent  deviations  from  the  seasonal  averages  of  these  earlier 
studies  and  may  fall  within  the  standard  deviations  of  the  seasonal  averages, 
which  were  not  discussed  by  the  earlier  authors. 


Significant  deviations  from  these  seasonal  averages  included; 

-  a  weak  undercurrent  in  the  presence  of  strong  equatorward  wind  stress  in 
May  1988; 

* 

-  the  lack  of  surface  equatorward  flow  during  a  period  of  strong  equatorward 
wind  stress  in  July  1989; 

-  the  presence  of  unusually  high  surface  salinities  of  southern  origin  in 

July  1989  suggestive  of  remote  forcing  from  the  south  resulting  in  unusually 
strong  poleward  flow  throughout  the  water  column; 

-  the  unusually  strong  and  deep  penetration  of  warm,  saline  waters  offshore  in 
November  1989  resulting  in  an  extremely  deep  signature  of  the  CC; 

-  generally  larger  magnitudes  observed  in  equatorward  flow  associated  with 
the  CC  resulting  from  a  deeper  LNM; 

-  larger  magnitudes  observed  in  poleward  flow  associated  with  the  CUC 
resulting  from  the  closer  spacing  of  stations  along  the  POST; 

-  definition  of  complex  banded  flow  near  the  coast  over  the  continental  shelf 
and  slope  resulting  from  the  closer  spacing  of  stations  along  the  POST; 

-  deep,  bottom  trapped  ageostrophic  flow  along  the  continental  slope  near 
2500  m  depth,  near  the  head  of  a  branch  of  the  Monterey  Canyon. 


This  last  feature  was  present  on  all  cruises  with  variable  intensity  due  to 
differences  in  station  spacing  and  the  variability  normally  associated  with  higher 
frequency  motions.  The  strength  and  persistence  of  the  phenomenon  suggest  that 
it  is  real  and  not  the  result  of  errors  in  sampling  and/or  computational  processes. 

B.  RECOMMENDATIONS 

This  study  focussed  on  only  seven  cruises  selected  to  represent  specific 
seasons.  A  better  understanding  of  seasonal  and  temporal  variability  would  be 
obtained  by  using  all  of  the  cruises  to  date  along  the  POST.  This  not  only 
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includes  using  the  hydrographic  data,  but  ADCP  and  PEGASUS  data  collected 

concurrently  with  CTD  data,  and  current  meters  which  have  been  placed  at 

several  locations  along  the  POST  since  1988.  The  time  scales  of  phenomena  off 
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Point  Sur  can  be  resolved  by  these  current  meters  and  will  aid  in  the 
interpretation  of  the  hydrographic  data.  As  this  data  base  continues  to  grow  so 
will  our  understanding  of  the  processes  which  occur  along  the  POST. 

The  wind  relaxation  and  subsequent  poleward  flow  in  July  1989  were 
suggestive  of  remote  forcing  from  farther  south.  Wind  records  from  NOAA 
buoys  to  the  south  were  not  available  at  the  time  of  this  study  however  they 
should  be  obtained  and  incorporated  into  future  studies  along  the  POST.  Such 
information  would  provide  insight  into  the  cause  of  such  intense  bursts  of 
poleward  flow  as  observed  in  July  1989. 

The  intrusion  of  warm,  salty  water  in  the  offshore  section  of  the  POST  was  a 
feature  common  to  several  cruises.  It  appears  to  be  either  North  Pacific  Central 
water  or  an  unusually  strong  recirculation  of  CUC  water.  A  more  rigorous  T-S 
analysis  is  required  to  determine  the  exact  cause  of  the  warm,  salty  water. 

Oceanography  data  frequently  requires  interpolation  to  provide  initial 
conditions  for  models,  estimation  of  transports,  and  wave  number  spectra 
(McIntosh  1990).  In  recent  years  the  use  of  objective  analysis  has  gained 
increasing  popularity  in  the  representation  of  hydrographic  data  (Bretherton  et 
al.  1976;  Roemmich  1983),  modeling  of  mesoscaie  features  (Reinecker  et  al. 
1987),  and  in  the  interpretation  and  analysis  of  satellite  imagery  (Kelley  and 
Caruso  1990;  Wahl  and  Simpson  1990).  Objective  analysis  commonly  refers  to  a 
statistical  interpolation  method  based  on  the  Gauss-Markov  theorem  (Bretherton 
et  al.  1976)  and  has  been  used  to  map  hydrographic  station  data  and  derived 
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fields  such  as  density,  salinity  anomaly,  and  geostrophic  velocity  to  standard 
positions  and  depths.  This  technique  provides  an  alternative  to  the  shear 
extrapolation  techniques  presently  used.  It  has  the  advantage  of  utilizing  all  the 
available  data  whereas  present  extrapolation  techniques  utilize  data  between 
adjacent  stations  pairs  to  fill  in  the  region  below  the  deepest  common  depth  to  a 
station  pair.  Such  techniques  are  artificial  and  may  result  in  larger  errors  than 
objective  techniques,  especially  when  computing  transports  (Roemmich  1983). 

Interpolation  of  temperature  and  salinity  data  to  a  uniform  grid  using  an 
objective  analysis  technique  similar  to  that  of  Bretherton  et  al.  (1976)  or 
McIntosh  (1990),  which  utilizes  splines  (a  curve,  properly  weighted,  so  that  it 
passes  smoothly  through  each  point)  would  provide  a  uniform  basis  upon  which 
more  rigorous  seasonal  and  interannual  comparisons  could  be  made.  This  would 
eliminate  problems  associated  with  stations  not  being  occupied  on  consecutive 
cruises  due  to  operational  considerations,  and  those  inherent  in  artificial 
extrapolation  techniques  where  the  bottom  depth  becomes  shoaler  than  the 
reference  level. 
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APPENDIX  A.  CALIBRATION  INFORMATION  AND  DEEP 

T-S  CURVES 

t 

The  following  scattergrams  and  regressions  were  the  result  of  the  calibration 
procedures  discussed  at  length  in  Chapter  III.  All  cruises,  with  the  exception  of 
those  discussed  earlier,  have  been  shown  here.  T-S  curves  for  the  2800  to  3300 
uecibar  range  have  also  been  included.  These  plots  consist  of  all  stations  along 
the  POST  that  had  depths  falling  into  this  range.  As  mentioned  in  Chapter  III, 
these  T-S  curves  were  used  throughout  the  calibration  process  to  determine  the 
effectiveness  of  the  calibrations. 
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August  1988  California  Undercurrent  Cruise  :  Large  Neil  Brown  CTD 
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November  1988  California  Undercurrent  Cruise  :  Small  Neil  Brown 
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118 


STMAY  1989  POST  Cruise  :  Large  Neil  Brown  CTD 
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True  Salinity  (PSU) 


July  1989  California  Undercurrent  Cruise  :  Large  Neil  Brown  CTD 
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APPENDIX  B:  PROPAGATION  OF  VARIANCES  AND 
COVARIANCES  AND  FUNDAMENTAL  PROPERTIES  OF  A 

UNIFORM  DISTRIBUTION 

The  material  presented  within  this  appendix  provides  a  detailed  description 
of  the  method  of  propagation  of  variance  or  propagation  of  errors  as  it  is  more 
commonly  referred  to,  following  Mikhail  (1976).  A  general  description  of  this 
process  will  be  provided  first,  followed  by  the  derivation  of  the  equations  found 
in  Chapter  IV. 

In  general,  observations  can  be  considered  samples  from  probability 
distributions  of  random  variables.  The  error  properties  of  these  observations 
are  then  considered  to  be  statistical  properties  of  the  sampling.  It  is  common  to 
assume  that  the  observations  are  statistically  independent.  By  this  it  is  meant  that 
the  results  of  a  previous  observation  does  not  or  should  not  influence  subsequent 
observations.  This  statistical  independence  results  in  the  elimination  of 
covariance  terms  between  random  variables.  If  the  observations  are  somehow 
correlated  with  each  other  they  are  no  longer  considered  statistically 
independent.  For  the  purposes  of  this  work  the  observations  are  assumed  to  be 
statistically  independent  both  horizontally  and  vertically. 

Propagation  of  error  involves  determining  the  random  characteristics  of 
dependent  variables  from  the  random  characteristics  of  the  independent  variables 
and  the  functional  relationships  which  relate  the  two  sets  of  variables.  Letting  xj 
and  X2  represent  random  variables,  |ii  and  p.2  the  means  and  f(\\,  X2)  the 
probability  density  function,  the  variances  and  covariances  may  be  defined  in  the 
following  manner: 
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(xi  - 


(*2  -  ^2?f{x2)  dx2 


B1 


B2 


crxix2  =  E[(xi  -  Hi  Xx2  -  |0-2 )] 

(Xi  -  PiX*2  -  Ji2if(xilf(x2)  dxidx2 

Bd 

If  y  l  —  gi(Xj,  x2)  and  y2  =  g2(xi,  x2),  then  the  following  expressions  for  the 
variances  and  covariances  for  yi  and  y2  are  obtained. 


(gt(xi,x2) -  Py^/fxi  ,x2)  dxjdx2 


(g2(xi,x2)-  Py^/Jxi  ,x2)  dx)dx2 


(gi(xi,x2)  -  pyi)(g2(xi,x2)  -  Py2If(x!  ,x2)  dxidx2 


B4 


B5 


B6 


General  formulae  such  as  these  are  seldom  used  in  practice.  The  prop¬ 
agation  of  variances  and  covariances  are  simplified  to  linear  or  linearized 
functions. 
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A.  PROPAGATION  OF  VARIANCE  AND  COVARIANCE  FOR 
LINEAR  FUNCTIONS 

Consider  the  following  linear  functions,  where  xi  and  X2  are  random 

* 

variables  with  means  |A]  and  p2  and  a  joint  probability  density  function  /(x i,X2), 


yi  =  ao  +  aiX!  +  a2x2 


B7 


y2  =  bo  +  bjXj  +  b2X2 


From  the  definition  of  variance  and  covariance  (B1-B3)  and  the  fact  that  both  the 
summation  and  integral  operators  are  linear,  expressions  for  the  variance  and 
covariance  of  yi  and  y2  are  obtained  in  the  following  manner. 

Oyi  =  E[(yi  *  Hy,)?] 

=  E[(a0  +  ajxi  +  a2x2  -  ao  -  aipi  - 
=  E[(ai(xi  -  pi)  +  a2(x2-p2))2] 

=  E[a?  (xi  -  pi  j2  +  (x2  -  P2)2  +  2aia2(xi  -  pi  \\2  -  P2)] 

=  ai  e[  (*1  '  Pi)2]  +  a2  E[  (x2  -  P2)2]  +  2aia2  E[(xj  -  Pi  Xx2  -  M-2>] 

=  afox,  +  a  22olz  +  2aia2  oX]X2  Bg 

Similarly  for  y2, 

Oy2  =  b](JXl  +  b2oX2  +  2b  1  b2  ox,x2  39 
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The  covariance  is  obtained  as  follows, 


Oyiya  =  Et(yi '  ^yi  Xy2  -  Ml 

=  E[(ao  +  aixi  +  a2X2  -  ao  -  aiPi  -  a2P2Xbo  +  bixi  +  b2X2  -  bo  -  biPi  -  b2ft2)] 

=  E[(a1(x1  -pi)  +a2(x2  -  ti2)Xbi(xi  -Pi)  +  b2(x2-p2))] 

=  Ejaib!  (xi  -  M-i)2  +  a2b2  (x2  -  P2 f  +  (aib2  +  aib2  Xxi  *  M-i  Xx2  -  P2 )  ] 

=  aibi  E[  (x!  -  (ii)2]  +  a2b2  E[  (x2  -  P2)2]  +  (aib2  +  aib2  )E[(xi  -  Pi  Xx2  - P2)] 

=  aibiO^  +  a2b2a^2 +  (aib2+ aib2  )aXlX2  BIO 


Upon  inspection  of  these  equations,  notice  that  they  are  independent  of  the 
density  functions.  Therefore  propagation  of  variances  and  covariances  of  linear 
functions  are  valid  for  any  distribution. 

If  matrix  notation  is  introduced  it  is  possible  to  consolidate  equations  B8-B10 
into  a  much  simpler  form.  Let  y  =  [yi  y2l 1  and  x  =  tx  1  x2l 1  »  where  x  and  y 
are  vectors  and  t  denotes  the  transpose  of  the  matrix.  Using  x  and  y  equation  B7 
may  be  rewritten  as, 

y  =  c  +  C  x  B 1 1 


where 


c  =  [  ao  bo  ]' 


C  = 


ai  a2 
bi  b2. 


B12 


4 
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Utilizing  equations  B8-B10  the  covariance  matrices  for  two  random  vectors 
x  and  y  can  be  written  as  follows, 

?  i  r  2 

CX1X2  _  °yi'  °yiy2 

^xx  -  ,  ^yy  ~  2 

_Oxlx2  °X2  J  l°yiy2  °y2  J  B13 

which  when  combined  with  B12  yields, 

£yy  =  C  SXx  C*  B14 

Here  the  matrix  C  represents  the  Jacobian  (J)  of  y  with  respect  to  x,  where  J  is 
written  as, 

J  =  ^ 

yx  3x  B15 


Finally,  substitution  of  B15  into  B14  the  general  form  of  the  propagation  of 
variance  and  covariance  becomes, 


£yy  —  jyx  ^xx  Jyx* 


B.  PROPAGATION  OF  VARIANCE  AND  COVARIANCE  FOR 
NONLINEAR  FUNCTIONS 

The  previous  discussion  dealt  with  linear  functions,  however,  in  practice  one 
is  usually  confronted  with  nonlinear  functions.  In  this  case  equations  such  as  B7 
are  the  result  of  a  linearization  of  the  general  equations.  If  we  now  consider 
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yi  =  gi(xi,  X2)  and  y2  =  g2(xi.  x2)  to  be  nonlinear,  the  linearized  form  at  the 
initial  values  xi°  and  X20  will  be  identical  to  B7. 

The  linearized  form  is  generally  truncated  after  the  first  order  terms,  thus 

r 

neglecting  any  higher  order  terms.  The  zero  order  terms  are  not  required  since 
gi(xi°,  X20)  and  g2(xi°,  X20  )  correspond  to  ao  and  bo  and  will  not  appear  in  the 
forms  for  propagated  variance  and  covariance  (B8-B10).  All  that  remains  are 
the  partial  derivatives  and  equation  B6  will  include  all  four  partial  derivative:. 


ayi 

dyi 

3xi 

3x2 

3v2 

3v'2 

3xi 

3x2  J 

evaluated  at  x  i°  and  X20. 

The  use  of  linearized  forms  is  common  practice  in  the  propagation  of 
variances  and  covariances  but  it  must  be  kept  in  mind  that  they  are  only  suited  to 
regions  where  the  function  is  well  approximated  by  its  tangent.  In  this  case,  the 
properties  of  the  nonlinear  random  variables  now  become  properties  of  the 
increments, 

Xi  =  xi0  +  Ax;  yi  =  yio  +  Ay, 

In  other  words,  the  error  properties  are  now  associated  with  Axj  and  Ayj  instead 
of  xi  and  yi,  respectively. 
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C.  *  ROPAGATION  OF  VARIANCE  THROUGH  THE  EQUATION 
OF  STATE  AND  EQUATION  FOR  GEOSTROPHIC  VELOCITY 

This  section  provides  the  details  involved  in  the  derivation  of  equation  5  in 

* 

Chapter  IV.  Following  the  methods  discussed  earlier  in  the  Appendix,  our 
function  becomes  5  =  S(S,T,P),  and  equations  B13  through  B16  become, 


^STP  = 


o|  Ost  Osp 
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From  statistical  independence  the  off  diagonal  terms  have  been  set  to  zero.  The 
final  result  becomes. 


B21 


which  is  equation  5  in  Chapter  IV. 

To  determine  propagated  errors  in  geostrophic  velocity  the  same  procedure 
shown  above  is  followed.  Starting  with  the  equation  for  geostrophic  velocity. 


(ADb  -  ADa)  +  vrcf 


B22 


9vg  dvg  dvg  dvg 

BADb  3ADa  ^Vrcf  .  B24 
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9ADb  3  AD  4 
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0  °ADa  0  0 

0  0  c?  0 
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9ADb 

avg 

9ADa 
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AD*  9L 


9vs,\2 


Substitution  of  the  derivatives  into  the  Jacobian  matrix  yields  the  following, 


«■ = (jtf  a^» + (tc)2  + (r  F  °l + fef  << 


B28 


Because  the  error  in  dynamic  height,  a^D2*  is  the  same  for  each  station,  the  first 
two  terms  in  equation  B28  are  equivalent  and  the  derivative  of  vg  with  respect  to 
vref  is  equal  to  one.  As  a  result  B28  reduces  to, 


(lFo£+<( 


B28 


This  form  is  identical  to  equation  (9)  in  Chapter  IV  where  the  value  in  (9) 
already  incorporates  the  factor  of  two. 

D.  UNIFORM  DISTRIBUTIONS 

This  section  contains  the  development  for  the  error  in  station  spacing  as 
described  in  Chapter  IV.  The  source  material  comes  from  Bendat  and  Piersol 
(1986). 

Consider  an  experiment  which  consists  of  choosing  a  point  at  random  in 
some  interval  [a,  b]  which  includes  the  endpoints.  The  probability  distribution 
for  a  continuous  random  variable  x(k),  becomes, 


p(x)  = 


’ 

0  x  <  a 

a  <  x  <  b  ; 

b  -  a 


1  x  >  b 
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and  the  corresponding  probability  density  function  is  defined  as. 


— - —  a  <  x  <  b 

b  -  a 

0  otherwise 


B30 
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The  mean  or  expected  value  is  defined  as, 


E[g(x(k))]  =  j  g(x)p(x)  dx 


B31 


and  for  this  example  becomes, 


b 


B32 


The  variance  or  second  moment  is  defined  as, 


E[x2(k)]  = 


x2p(x)  dx 


B33 


and  for  this  example  becomes. 


o 


X 


(b-a)2 

12 
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For  the  case  of  the  error  in  station  spacing  due  to  vessel  drift,  the  value  of 
the  error  is  assumed  to  lie  between  0  and  L/2,  or  V2  the  differential  drift 
between  stations.  Substituting  these  values  in  for  a  and  b  in  equations  B19 
through  B22,  the  following  result  is  obtained. 


B35 


Taking  the  square-root  of  this  result  provides  the  root  mean  square  error  which 
for  the  error  in  station  spacing  is  simply  the  differential  distance,  L/2,  divided  by 
two  times  the  square-root  of  three. 
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